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ABSTRACT 


•  1 

-  Optical  variability  across  the  continental  slope  and 
shelf  off  Central  California  was  stalled  asing  Nimbus-7 
Coastal  Zone  Color  Scanner  (C2CS)  data.  CZCS  estimates  of 
k(490),  the  irradiance  attenuation  coefficient  at  490  nm, 
were  expressed  as  optical  depth  1/k(49Q).  A  modified  atmos¬ 
pheric  correction  algorithm  was  used  to  account  for  water 
radiance  at  670  nm.  Time  sequences  of  1/k(490)  were  assem¬ 
bled  and  partitioned  into  four  zonal  transects,  at  different 
latitudes,  spanning  Hay  through  November  in  1979,  1980  and 

1982.  Empirical  Orthogonal  Functions  (EOFs)  were  calculated 
for  each  partition.  The  first  EOFs  are  dominated  by  scales 
of  order  180  km,  with  in  all  cases,  a  band  of  low  optical 
depth  water  in  the  first  100  km  adjacent  to  the  coast. 
Scales  decrease  in  successive  EOFs,  to  about  40  km  in  the 
fifth  EOF.  The  feasibility  of  joining  EOFs  from  different 
partitions  was  demonstrated  as  a  precursor  for.  future 
applications  to  piecewise  analysis  of  oceanic  satellite 
da  ta . 

/V 
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I.  INTRODUCTION 

Satellite  remote  sensing  systems  offer  fast,  eccr.cnicai 
means  of  determining  the  horizontal  structure  of  the  oceans 
on  a  global  basis.  The  objective  of  this  thesis  is  to  con¬ 
tribute  to  the  development  of  empirical  methods  for  using 
satellite  images  of  optical  parameters  and  sea  surface  temp¬ 
erature  (SST)  to  infer  the  upper  ocean's  vertical  structure, 
through  interpolation  and  extrapolation  of  relatively  lim¬ 
ited  in  situ  data. 

The  objective  is  being  approached  through  regional  case 
studies  of  correlations  between  optical  parameters  and  phys¬ 
ical  water  mass  properties  in  the  upper  ocean  in  different 
regions  of  the  world.  Mere  specifically,  this  thesis  is  a 
preliminary  case  study  of  the  California  Current  region. 
The  ultimate  goal  is  to  relate  statistically  the  horizontal 
structure  cf  optical  properties  observed  with  the  Coastal 
Zone  Color  Scanner  (CZCS)  to  the  underlying  vertical  struc¬ 
tures  of  temperature  and  salinity,  as  well  as  bio-optical 
parameters,  for  a  given  region  and  season. 

The  study  domain  encompasses  the  continental  slope  and 
shelf  off  the  coast  of  California  between  Point  Sur  and 
Point  Arguello.  This  area  was  selected  to  investigate  an 
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ocean  upweixmg  rror.t  vnrcn  ns  kr.cvn  to  persist  throughout 
the  upwallir.g  season  (Traganza,et  al. ,  1979).  The  northern 
and  southern  portions  of  the  study  domain  are  typified  by 
complex  eddy  structure  associated  with  irregular  features  in 
the  bathymetry,  such  as  off  Point  Sur.  Between  Point  Sur 
and  Point  Arguellc,  on  the  other  hand,  isolines  of  S3T  and 
optical  parameters  tend  to  be  aligned  roughly  parallel  to 
the  underlying  isobaths. 


An  ensemble  of  data  acquired  with  the  Nimbus  7  CZCS  dur¬ 


ing  the  summer  and  fall  seasons  of  1979,  1980  and  1982  is 
analyzed  in  this  study.  The  horizontal  structure  in  bio-op- 
tical  parameters  determined  from  cloud-free  portions  of  CZCS 


imagery  are  investigated  using  a  Partitioned  Empirical 
Orthogonal  Function  (PEOF)  decomposition.  The  spatial  par¬ 
titions  examined  here  consist  of  four  zonal  transects  cross¬ 


ing  the  shelf/slope  region  at  different  latitudes.  The 
specific  goals  of  this  analysis  are: 

1.  To  characterize  the  meridional  and  zonal  spatial 
correlation  structures  of  ocean  color  parameters 
(specifically  optical  depth  1/K(490)  in  meters). 


2.  Tc  compare  the  spatial  scales  and  structures  of 
optical  variability  highlighted  by  the  PEOF 


5 


A.  THE  STUDY  DOMAIM 
1 .  Coverage 

a.  Area  Domain 

The  region  investigated  in  this  project  is 
located  between  32  and  40N,  and  from  the  coast  of  California 
offshore  to  approximately  126W,  Pig.  1.  This  area  was 
selected  because  it  contains  water  mass  structures,  includ¬ 
ing  fronts,  which  strongly  influence  phytoplankton  concen¬ 
trations,  and  therefore  the  optical  properties  of  the  ocean 
water.  Furthermore,  an  adequate  sample  of  data  was  availa¬ 
ble  for  this  area. 

A  subarea  of  this  region  is  labelled  Insert  A  in 
Figure  1  and  presented  in  greater  detail  in  Figure  2  Insert 
A  is  bounded  by  34  to  38N,  and  by  126  to  120w.  It  is  the 
primary  study  domain  of  this  thesis.  The  background  hydrog¬ 
raphy  and  dynamics  of  this  region  are  described  in  subsec- 


b.  Time  Domain 

In  the  time  domain,  the  available  CZCS  data 
include  scenes  from  summer  through  early  fall  seasons  in 
1979,  1980  and  1982.  Originally,  a  single  season  ensemble 
of  CZCS  data  (May  through  September  1980)  was  sought.  How¬ 
ever,  the  time  span  had  to  be  expanded  to  three  years  to 
assemble  a  reasonably  large  sample  size  of  cloud-free  sub¬ 
scenes.  The  sample  analyzed  contains  three  scenes  in  1979, 
eleven  scenes  in  1980  and  eight  scenes  in  1982.  Detailed 
characteristics  of  these  CZCS  images  are  presented  in 
Appendix  B. 

2 .  Geometry 

a .  Coastal 

The  California  coast  bounding  the  area  of  study 
is  characterized  by  a  steep,  mountainous  coastal  range  run¬ 
ning  roughly  parallel  to  the  coastline.  The  coastline 
stretching  from  San  Francisco  to  Point  Arguello  is  oriented 
roughly  northwest  to  sout least,  but  is  interrupted  by  Monte¬ 
rey  Bay  at  36  45'N  and  by  smaller  bays  in  the  vicinity  of 
Morro  Bay  at  35  20  'N.  No  major  rivers  drain  into  this 

coastline,  although  many  local  rain-generated  drainage 


creeks  empty  here. 


Bathymetry 


b . 

The  predominant  orientation  of  the  bathymetry  is 
roughly  northwest  to  southeast  (parallel  to  the  coast).  Fig. 
2.  Interruptions  of  this  orientation  are  evident  in  the 
vicinity  of  the  Monterey  Canyon,  Point  Sur  and  the  Sur  Can¬ 
yon,  the  Davidson  Seamount,  the  Taney  Seamount  and  the  Santa 
Lucia  Banks  and  Escarpment  off  Point  Arguello.  k  very 
abrupt  shelf  break  is  evident  all  along  this  section  of  the 
California  coast.  Isobaths  tend  to  diverge  south  of  Monte¬ 
rey,  due  to  a  broadening  of  the  continental  shelf  and  slope 
with  distance  south  of  Monterey. 

3.  Descriptive  Oceanography 
a.  Coastal  Opwelling 

Coastal  upwelling  is  an  oceanic  phenomenon  which 
has  a  pronounced  impact  upon  many  physical  and  biological 
processes.  Predominantly  southward  winds  during  spring  and 
summer  off  the  central  California  coast,  yield  offshore  sur¬ 
face  Ekman  transports,  which  forces  compensation  water  to 
rise  from  depths  of  the  order  of  200  to  300  m  (Smith,  1968). 

The  upwelling  season  off  the  coast  of  California 
is  generally  confined  to  the  late  spring  through  early  fall. 
The  onset  of  the  seasonal  upwelling  commences  in  more 
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southern  waters  off  the  California  coast  and  progresses 
northward  as  the  season  unfolds  (Yoshida  and  Mao,  1957; 
Wooster  and  Reid,  1963;  Pavlova,  1966;  Hickey,  1979).  The 
CZCS  data  set  spans  the  upwelling  season  and  includes  images 
from  beyond  this  season  into  early  winter. 

Opwelling  has  a  marked  effect  on  the  sea  surface 
temperature,  causing  it  to  be  much  lower  than  would  other¬ 
wise  be  normal  for  the  latitude  and  season  (Smith,  1968). 
The  relatively  lower  temperatures  are  evident  in  IR  images 
of  the  region  (Johnson,  1980;  Nestor,  1979).  Accompanying 
this  decrease  in  temperature  is  an  increase  in  surface 
salinity,  an  upwelling  property  unique  to  the  regime  off  the 
west  coast  of  North  America  (Smith,  1963)  . 

The  oceanographic  properties  of  upwelling  have 
been  documented  in  many  areas  of  the  world,  but  nowhere  with 
the  thoroughness  of  the  work  off  California  and  Oregon. 
Ship  and  satellite  observations  have  allowed  us  to  identify 
seasons,  centers,  and  the  extent  of  the  upwelling  event 
along  the  west  coast  of  Ncrth  America.  Traganza,  et  al. 
(1  979)  used  combined  satellite  and  shipboard  observations  to 
infer  nutrient  upwelling  distributions  off  the  coast  of  Cal¬ 
ifornia.  Frontal  structures  and  mesoscale  eddies  that  can 


result  from  the  upwelling  phenomenon  have  been  examine  d  with 
relevance  to  Anti-Submarine  Warfare  (ASW)  by  Tragenza 
(1979).  The  use  of  infrared  (IE)  imagery  in  the  detection 
and  description  of  upwelling  was  examined  by  both  Johnson 
(1980)  and  Nestor  (1979). 

The  introduction  of  nutrient  rich  waters  to  the 
nearshore  euphotic  zone  greatly  enhances  the  development  of 
the  in  situ  phytoplankton  population.  This  enhancement  in 
turn  causes  the  upwelled  water  mass  moving  offshore  to  have 
distinctly  different  optical  properties  than  adjoining  off¬ 
shore  waters.  The  boundary  (frontal  region)  between  the 
upwelled  water  mass  and  the  normal  surface  water  mass  is 
thus  readily  detectable  and  of  great  interest. 

Nutrient  enrichment  off  the  California  coastal 
zone  is  observed  in  the  regions  of  upwelling  evenns.  These 
nutrients,  which  are  classified  as  "biochemically  new”  on 
the  basis  of  nitrat e-t o-phosph ate  ratios  which  approach 
15:1,  are  brought  to  the  surface  from  depths  up  to  300  m. 
By  way  of  contrast,  nutrients  also  present  in  the  open  ocean 
surface  water  approach  5:1  (Nestor,  1979).  The  added 
nitrates  are  a  primary  factor  in  the  increase  in  phytoplank¬ 


ton  concentrations  during  the  upwelling  season 


Coastal 
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opa?.  ocean  waters  and  this  difference  can  be  increased  dur¬ 
ing  periods  of  upwelling  (Sverdrup,  et  al.,  1  942).  The  phy¬ 
toplankton  concentrations,  with  their  associated 
chlorophyll-like  pigments,  have  a  profound  effect  or.  +he 
upwelled  radiances  measured  by  the  CZCS,  as  discussed  in 
Chapter  III. 

Another  aspect  of  upwelling,  and  its  relation  to 
satellite  data,  is  its  effect  on  regional  climate.  The  rel¬ 
atively  cold  sea  surface  temperature  in  upwelling  zones 
cools  the  air  above  and  thus  increases  its  relative 
humidity.  As  a  result,  low  stratus  and  fog  commonly  occur 

hare  in  a  shallow  (marine)  layer  with  warm  air  aloft.  The 
frequent  occurrence  of  low  stratus  and  fog,  seriously  limits 
infrared  and  visible  satellite  coverage  during  the  upwelling 
season.  The  cool  sea  water  also  contributes  to  a  diurnal 
sea  breeze  by  increasing  the  onshore-offshore  pressure  gra¬ 
dient.  Onshore  wind's  bring  cool,  moist  air  as  far  as  50 
miles  inland  (Smith,  1968)  . 
b.  Currents 

The  California  Current  System  may  be  discussed 


and  studied  in  terms  of  four  large  scale  currents: 


the  California  Current,  the  California  Undercurrent,  rh? 
Davidson  Inshore  Current,  and  the  Southern  California  Cur¬ 
rent  (Hickey,  1979).  The  first  three  of  these  currents 
directly  influence  the  study  domain.  Masoscale  currents 
associated  with  seasonal  upwelling  are  also  important  here. 

(1)  California  Cur£ er^t.  The  California  Current 
is  a  broad  wind-driven  equatorvard  current  which  exhibits 
significant  seasonal  variations  proportional  to  the  changes 
in  the  wind  field  (Brown,  1974)  .  Off  Point  Conception  the 
mean  annual  location  of  the  current  axis  is  located  270  km 
offshore  while  the  shoreward  boundary  extends  to  200  km  off¬ 
shore.  The  current  is  cf  the  order  700  km  wide  and  flows 
south  at  10  to  30  cm  per  second  (Hickey,  1979). 

The  Califcrnia  Current  is  a  continuation  of 
the  West  Wind  Drift  in  the  North  Pacific  and  flows  southward 
along  the  California  coast  between  48  and  23N.  It  turns 
westward  between  20  and  30N  where  it  becomes  part  of  the 
North  Equatorial  Current.  This  flow  regime  comprises  the 
eastern  extent  of  the  anticyclonic  HE  Pacific  Subtropical 
Gyre,  which  is  centered  near  the  Hawaiian  Islands  (Sverdrup, 
et  al.,  1942;  Chelton  and  Davis,  1982) . 


(2)  California  Under current.  Ine  California 
Undercurrent,  also,  referred  tc  as  the  California  Ccun-er- 
current,  is  the  poleward  subsurface  flow  over  the  continen¬ 
tal  slope.  Maximum  poleward  flew  occurs  during  the  summer 
and  fall  seasons  at  depths  of  200  to  250  m  (Pavlova,  1966 
and  Hickey,  1979).  The  flow  can  be  described  as  a  broad 
current  with  a  central  jet.  It  is  this  jet  structure  that 
is  most  eften  measured  and  referred  to  when  applying 

specific  values  to  the  Undercurrent.  The  broad  poleward 

flow  has  a  gecstrophic  component  alongshore  near  the  shelf 
break  of  approximately  15  cm/sec  (Ccddington,  1979). 

The  flow  appears  to  have  a  jet-like  struc¬ 
ture,  both  vertically  and  horizontally,  and  to  extend  to  the 
bottom  over  the  slope.  The  existence  of  a  high  speed  jet 
core  of  the  order  of  20  to  70  km  in  width,  was  first  sug¬ 
gested  by  Beid  (1962,  1963).  Subsequent  direct  measurements 
of  these  jets  have  produced  values  as  high  as  40  cm/sec  off 
Northern  Baja  and  values  of  16  cm/sec  off  Washington 
(Wooster  and  Jones,  1  970)  .  The  depth  of  the  high-speed  core 
varies  seasonally.  It  rises  from  depths  of  200  to  300  m  to 

the  surface  during  the  late  fall  and  winter  north  of  Point 


Conception 


Here  it  is  referred  as  the  Davidson  Inshore 
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Current  by  many  authors  (Hickey,  1979;  Pavlova,  1966;  Ingra¬ 
ham,  1967).  Event-scale  fluctuations  (of  rhe  order  of  ICO 
km  and  10  days)  in  the  flew  appear  to  be  correlated  with  the 
alongshore  component  of  wind  stress  (Nelson,  1977)  . 

The  extent  and  time  scale  of  continuous 
alongshore  flow,  and  the  width  of  the  region  of  northward 
flow  below  500  m,  are  important  topics  yet  to  be  answered 
about  the  California  Undercurrent  (Hickey,  1979). 

(3)  Davidson  Inshore  Current.  North  of  Point 
Conception,  the  poleward  surface  flow  in  the  nearshore 
ragiens  off  the  West  Coast  is  known  as  the  Davidson  Inshore 
Current.  It  is  associated  with  winter  weather  circulation 
patterns.  As  the  southward  winds  weaken  and  tend  toward  a 
northwestward  flow,  the  Davidson  Inshore  Current  becomes 
established  (Hickey,  1979)  .  The  current  flows  near  the 
coast,  usually  within  100  km,  well  inshore  of  the  California 
Current  and  is  confined  tc  the  continental  shelf  and  slope. 
Pavlova  (1966)  reported  that  north  of  Point  Conception,  the 
Davidson  Inshore  Current  reaches  its  maximum  development  at 
depth  (200  to  250  m)  in  the  summer  and  autumn.  In  August 
the  Davidson  Inshore  Current  is  scarcely  noticeable  at 
the  surface  despite  active  development  a  depth.  Maximum 
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surface  development  is  reached  from  October  through  Ap ril, 
i.  a.  ,  late  autumn  to  early  spring.  In  December,  the  cere  of 
maximum  velocity  emerges  at  the  surface  and  in  the  late 
spring  it  almost  completely  disappears  (Reid,  1960;  Reid, 
et  al.,  1958;  Favlova,  1S66).  Poleward  velocities  of  up  to 
25  cm/sec  were  recorded  (Eeid  and  Swartzlose,  1962)  within 
80  km  of  central  California  in  January. 

The  Davidson  Inshore  Current  and  the  Cali¬ 
fornia  Undercurrent  are  often  discussed  as  though  they  were 
separate  currents.  Both  currents  transport  Equatorial-type 
water  northward  at  least  as  far  as  Cape  Mendicino  (Pavlova, 
1966).  Also,  no  subsurface  maximum  has  been  found  in  the 
flow  of  the  Davidson  Inshore  Current.  These  characteristics 
support  a  view  that  the  Davidson  Current  is  simply  the  sur¬ 
face  expression  of  the  California  Undercurrent,  rather  than 
a  separate  current  superimposed  cn  it. 

(4)  Other  Currents.  The  presence  of  eddies 
throughout  the  the  California  Current  System  has  been  docu¬ 
mented  for  many  years  (Bernstein,  et  al.  ,  1977).  The  time 

scales  for  these  eddies,  as  well  as  the  processes  responsi¬ 
ble  for  their  generation  and  subsequent  dissipation,  is  an 
area  of  increasing  study. 
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Between  San  Francisco  and  a  point  about 
half  the  distance  to  Point  Conception,  there  is  a  permanent 
counterclockwise  eddy  that  produces  northward  flow  during 
all  months  except  April  (Erown,  1974).  A  second  eddy  just 
north  of  Point  Conception  forms  during  the  summer  months  and 
makes  northward  flow  continuous  from  Point  Conception  to  San 
Francisco  (Hickey,  1  979).  willmott  (1983)  has  shown  that 
these  features  may  be  produced  by  flow  separation  of  the 
California  current  in  the  vicinity  of  major  coastal  capes. 

Reid,  at  al.  (1963)  made  direct  measurements  of 
an  eddy  (90  km  in  diameter)  off  the  northern  coast  cf  Baja 
California.  Hypotheses  for  eddy  formation  discussed  in 
their  paper  are  as  follows: 

(1)  The  process  of  upwelling  and  the  offshore  movement  of 
the  cclder,  more  saline  waters  might  degenerate  into  eddies. 
The  lateral  shear  between  the  upwelling  flow  away  from  the 
coast  and  the  California  Current  and  Undercurrent  (baro- 
tropic  instability)  could  produce  eddy  structures.  Tempera¬ 
ture  and  salinity  differences  set  up  strong  baroclinic  zones 
along  the  upwelling  boundaries  which  could  result  in  eddy 
formation.  (Sverdrup  and  Fleming,  1941) 
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(2)  The  offshore  surface  flow  during  upvelling  should  pro¬ 
duce  a  counter  current  (Munk,  1950) •  If  there  is 
substantial  north-south  variation  in  the  intensify  of  the 
winds,  then  seperate  countercurrents  of  different  strengths 
might  occur  along  the  coast. 

(3)  A  second  hypothesis  proposed  by  Reid,  et  al.  (1963),  is 
that  the  deeper  countercurrent  may  transfer  momentum  upward 
to  the  surface  layers,  at  times  when,  cr  in  regions  where,  a 
surface  current  does  not  prevail.  This  could  cause  spot 
intrusions  of  colder  circulating  waters  that  form  eddies 
where  neither  surface  countercurrents  nor  coastal  upwellirg 
produce  them. 

Additionally,  the  effects  of  bottom  steering  by  coastal 
topography,  and  the  associated  trapped  motions  must  be  con¬ 
sidered  when  discussing  eddy  formation.  (Hurlburt,  1974; 
Johnson,  1982;  Willmott,  1983) 
c.  Water  Hasses 

Descriptions  of  the  water  masses  that  contribute 
to  the  California  Current  System  are  given  in  Tibby  (1941), 
Sverdrup,  et  al.,  (  1942)  and  Reid,  et  al.,  (1958).  Four 
major  sources  are  discussed  by  the  authors; 

(1)  Subarctic  Water  Mass  -  from  the  north. 

(2)  Central  Water  Mass  -  from  the  west  and  northwest. 


(3)  Equatorial  Water  Mass  -  from  the  south. 

(4)  Water  derived  from  up  welling  sources. 

These  sources  were  simplified  in  Tibby  (1941)  and  Sverdrup, 
et  al.,  (1942)  into  two  extreme  sources  named  "Subarctic 
North  Pacific"  and  "Equatorial  Pacific". 

The  percent  of  each  water  mass  comprising  a  sam¬ 
ple  can  be  defined  by  entering  Figure  3  with  a  T-S  pair. 
However,  the  determination  of  percentage  composition  by  this 
means  cannot  be  used  for  water  above  depths  of  about  100  m. 
This  restriction  is  due  to  vertical  mixing  in  the  nearsur¬ 
face  layer  related  to  the  effects  of  wind  and  local  changes 
due  to  heat  and  mass  fluxes  across  the  air-sea  interface. 
Any  mixing  along  surfaces  of  constant  at  would  be  severly 
masked  in  these  shallower  depths  by  the  effects  of  turbulent 
vertical  mixing.  Also,  below  1000  m  the  differences  in  the 
T-S  relationships  of  the  two  extreme  water  masses  are  negli¬ 
gible.  For  intermediate  depths,  as  might  be  expected,  the 
percentage  of  equatorial  water  decreases  in  the  direction  of 
northward  flow.  The  Undercurrent  is  characteristically 
warmer  and  more  saline  than  the  California  Current,  and  it 


has  a  salinity  maximum  cn  the  a 


26.54  surface.  Off 


III.  CZCS  OCEAN  COLOR  IMAGES  AMD  UPPER  OCEAN 

- g?Tirar~Ptro?TFTir5 - 

A.  INTRODUCTION 

The  physical  processes  of  absorption  and  scattering 
relate  the  upweliing  radiance  just  beneath  the  sea  surface 
to  the  constituents  of  the  water  (Gordon,  1976)  .  Except  for 
coastal  waters  and  waters  influenced  by  river  discharge, 
biological  constituents  play  a  dominant  role  in  these  pro¬ 
cesses  (Smith  and  Baker,  1978;  Jerlov,  1976).  Optically,  the 
most  important  biological  constituent  is  phytoplankton, 
microscopic  plant  organisms  that  photosynthesize  and  make  up 
the  first  link  of  the  oceanic  food  web  (Steele,  1970). 
Ch lorophyll-a  is  the  dominant  photosynthetic  pigment,  and 
absorbs  light  strongly  in  the  blue  and  red  regions  of  the 
visible  spectrum  (400  to  700  nm)  (Hovis,  et  al.,  1980). 
Therefore,  as  the  concentration  of  phytoplankton  increases, 
the  color  cf  the  water  is  shifted  toward  green  hues  from  the 
deep  blue  cf  its  pure  state.  By  measuring  upwelled  radiance 
(backscattered  daylight)  in  specific  spectral  bands,  we  can 
determine  the  concentrations  of  phytoplankton  pigments  in 
the  ocean  (Gordon,  et  al.  ,  1980;  Gordon,  et  al. ,  1983). 


This  chapter  first  describes  the  CZCS  sensor  ar.  i  i-.s 
capabilities,  and  then  the  measured  signal  is  discus- ed. 
Algorithms  that  are  currnetly  applied  to  this  signal  to  cor¬ 
rect  for  atmospheric  effects  are  discussed.  Finally,  the 
algorithms  designed  to  convert  the  corrected  radiance  values 
to  phytoplankton  concentrations,  C,  and  irradiance  attenua¬ 
tion  coefficient,  k,  are  presented. 

3.  SYSTEM  DESCRIPTION 

1.  The  Nimbus- 7  Coastal  Zoqe  Color  Scanner  (CZCS) 

The  CZCS  was  built  by  the  Ball  Brothers  Research 
Corporation  tc  NASA's  specifications.  The  instrument  is  a 
spatially  imaging  multispectral  scanner.  Six  spectral  bands 
are  precisely  coregistered  and  internally  calibrated.  The 
swath  width  cf  the  CZCS  is  slightly  more  than  1600  km. 
Characteristics  of  its  five  visible  (443,  520,  550,  670,  750 
nm)  and  one  thermal  13  (10.5  to  12.5  ym)  channels  are  summa¬ 
rized  in  Table  I.  The  CZCS  has  an  active  scan  of  78  degrees 
centered  cn  nadir  and  a  field  of  view  of  0.0485  degrees, 
yielding  a  geometric  instantaneous  field  of  view  of  825  m 
(at  nadir)  from  a  spacecraft  altitude  of  955  km.  It  can 
tilt  the  scan  plane  20  degrees  from  nadir  in  2  degree  incre¬ 


ments  along  the  satellite  track  to  minimize  the  influence  of 
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su: 


direct  sun  glint. 


The  Nimbus- 7  spacecraft  is 


synchronous  orbit 


with  ascending  node  near 


local  r.cor. . 


TABLE  I 

Character istics  of  the  CZCS 


(Ho v is,  et  al.,  1980) 


Band 

Number  Wavelength 
<nm)  * 


1 

433 

to 

453 

2 

510 

to 

530 

3 

540  to 

56  0 

4 

660 

to 

680 

5 

700 

to 

800 

6 

10,500 

to 

12,  500 

Gain 


1 

0 

3 

2 

1 

0 


1 

0 

3 

2 

1 

0 


Saturation 

Measured 

Radiance 
(mW/cm2sr  pm) 

signal/noise 

5.41 

7.64 

9.23 

11.46 

158/1 

3.50 

5.10 

6.20 

7.64 

200/1 

2.86 

4.14 

5.10 

6.21 

176/1 

1.34 

1.91 

2.32 

2.88 

11&/1 

23.90 

350/1 

0.22  K* 

*  Noise  equivalent  temperature  difference  at  270  K. 


scattering  by  the  air  (Rayleigh  scattering)  and  by  micro 

scopic  particles  suspended  in  the  air  (aerosol  scattering) 

both  cf  which  increase  the  radiance  detected  at  the  sensor. 

Fresnel  reflecticn  (sun  glint)  can  be  ignored  as  the  tilting 

capability  of  the  CZCS  minimizes  its  effect.  However  the 

scattering  effects,  both  Rayleigh  and  aerosol,  must  be 

removed  from  I  (X)  to  give  usable  values  for  the  upwelled 

radiance  L  (X). 

w'  ' 

C.  CZCS  GEOPHYSICAL  ALGORITHMS 
1 .  Atmospheric  Corrections 

Frcm  the  signal  description  in  section  A  of  this 
chapter,  we  can  construct  the  following  formula 


Lt(X)  =  Lr(X)  +  La(X)  +  td(X)Lw(X) 

where 

L  t  =  Total  radiance 

L  r  =  Radiance  due  to  Rayleigh  scattering 
La=  Radiance  due  to  aerosol  scattering 


(1) 


Lw=  Upwelled  radiance  from  beneath  the  sea  surface 


td=  Diffuse  transmittance  of  the  atmosphere 


X  =  Wavelength 
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As  previously  mentioned  Lt  is  the  total  radiance 
measured  fcy  the  CZC3.  The  Rayleigh  scattering  tern  car.  be 


sxpressel  as 

F  (A)  t  (A) 


L  (A)  = 
r 


where 


u;  r 

- -  P  (a 

„  r 
o  s  9  L 


)  +  {  P  ( 0 )  +  P  ( 0  )  )P„( a  )  T  (A)  , 


- 


0  3 


(2) 


F  =  The  instantaneous  extraterrestial  solar 

o 

irradiance. 


tr  *  The  Rayleigh  optical  thickness  of  the  atmosphere. 

Pr  =  The  Rayleigh  scattering  phase  function. 

o_  =  The  scattering  angle  through  which  photons  are 
backscattered  from  the  atmosphere  to  the  sensor 
without  interacting  with  the  sea  surface. 

a+  =  The  forward  scattering  angle  of  those  photons 
which  are  scattered  in  the  atmosphere  toward  the  sea 
surface  (sky  radiance)  and  then  specularly  reflected 
frcm  the  surface  into  the  field  of  view  of  the  sen¬ 
sor  (p(9  )  term)  as  well  as  photons  which  are  first 
specularly  reflected  from  the  sea  surface  and  then 
scattered  by  the  atmosphere  into  the  field  of  view 
of  the  sensor  ( p  ( Q )  term). 


p 


The  Fresnel  reflectance  cf  the  air-sea  interface 


T  =  The  two-way  ozone  transmittance  of  the  atmosphere. 

03 

e  =  The  sensor  zentith  angle  at  me  observed  point  or. 
the  sea  surface. 

0  =  The  solar  zenith  angle  at  the  observed  point  on 

0 

the  sea  surface. 

The  aerosol  scattering  term  is  found  using  the  670 
nm  channel,  where  there  is  only  a  negligibly  small  contribu¬ 
tion  by  the  Lw  term.  (This  is  referred  to  as  the  "black 
ocean"  assumption.)  We  calculate  aerosol  radiance  at  X  * 
67  0  as 

L  (670)  =  L  ( 6  7  0 )  -  L  (670)  .  (3) 

ci  w  r 

The  key  assumption  in  this  algorithm  is  that  the  ratio  of 
aerosol  wavelengths  is  constant  over  a  scene,  and  is  given 
as 

L  (A) 

e  (  A  ,  6  7  0 )  =  — -  *  (4) 

L  (670) 
a 

e  is  calculated  using  either  simultaneous  direct  radiance 
measurement  from  a  ship,  or  upwelled  radiance  values  mod¬ 
elled  at  a  clear  water  pixel  (Gordon  and  Clark,  1981)  .  The 
latter  method  is  discussed  in  section  2  of  this  chapter. 
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Returning  to  equation  (  1)  , 


(A)  is  the  diffuse 


d 

'Transmit*: ance  of  the  atmosphere  ana  sea  surface,  which  may 
be  approximated  as 
[1  -  p (0)  I 

td(A)  =  - '£ -  exp 

m 

where  all  terms  have  been  previously  defined  except  m ,  which 
is  the  index  of  refraction  of  water  relative  to  air  and  is 
assumed  to  be  4/3  for  the  wavelengths  (400  -  700  nm)  . 

We  have  now  developed  the  basics  for  extracting  the 
upwelled  radiance  values,  from  the  CZCS  detected  signal, 

Lf 

2.  Clear  Water  Radiance 

The  scene  constant, e  ,  given  in  equation  (4)  is  cal¬ 
culated  as 


[  T  r  (  A  ) 

— -  +  T  (A) 

2  °3 

/  cos  9 

F  (A)  exp[-x.  (A)(sec9  +  sec0  )] 

e(A,A  )  -  U/A  >n  -2 - ^ - 2 - 

F  (A  )  exp[-Tn  (A  ) (secS  +  sec0  )] 
o  o  0  3  o  o 


where  r.  is  called  the  Angstrom  coefficient.  Equation  (6) 
can  te  rewritten  as 


F  ( A )  T  (A) 
n  0  0  3 

£ ( A , A  )  =  (A/A  ) 

0  0 


F  ( A  )  T  n  (A  ) 


(7) 


I 

N 

j 


i 


i 


where  T  n  (X)  =  €xpj^-TQ_(X  )  (sec  9  +  sec9o  )J  and  where  for 

the  CZCS  Xq=  670  nm. 

Gordon  and  Claric  (1981)  developed  the  concept  of 
clear  water  radiance  for  atmospheric  correction  of  CZCS 
imagery.  The  strategy  employed  in  this  study  was  to  find  an 

area  cf  the  image  that  could  be  assumed  -o  have  a  chloro- 

3 

phyll  concentration  less  than  0.25  mg/m  .  it  this  low  con¬ 
centration,  L  at  520  and  550  nm  are  assumed  to  be 
essentially  constant  for  a  given  solar  elevation.  Then, 
given  these  "clear  water  values"  cf  LW(A)  at  one  position, 

L  (x)  is  calculated  using  equation  (1).  e(X,670)  are  found 

a 

from  equations  (4)  using  the  computed  L  a(  X)  and  La(670) 
value.  Finally,  rearranging  equation  (7)  we  find  that 


In 

e  (  A  ,  A  )  / 

f  F  (A)  Tn  (X)  1 

o  O  3 

j 

n  (  X )  = 

0 

F  ( A  )  T  (A  ) 

m  O  O  U  3  0  _ 

i 

In  ( X / X  ) 

0 

Values  for  n  at  520  and  550  nm  (n(520)  and  n(55Q))  are  com¬ 
puted,  then  averaged  to  estimate  n(443).  The  Angstrom  coef¬ 


ficient  at  443  nm  cannot  be  directly  measured  in  this  way 


because  (442)  is  highly  sensitive  to  aven  minute  fluctua¬ 
tions  in  chlorophyll  concentrations  at  low  concentration. 

An  important  aspect  of  this  algorithm  is  that  nei¬ 
ther  surface  measurements  of  Lw  (A)  ,  nor  any  properties  of 
the  aerosol  are  reguired  to  implement  the  atmospheric 
correction. 

3 .  Bic^optic  Parameters 

a.  Chlorophyl  Concentrations 

Determination  of  chlorophyll  concentra tions 

C  from  ratios  of  L  ( A )  relates  the  surface  value  of  C  to  the 

w 

ratio  of  the  upwelled  radiance  at  two  different  wavelengths 
(Morel  and  Prieur,  1977;  Gordon  and  Clark,  1981).  The  basis 
for  this  is  that  to  a  first  approximation  Lw  is  proportional 
to  the  ratio  of  the  volume  backscattering  coefficient, 

B  (A  )  b  (A)  ,  and  the  volume  absorption  coefficient,  a  (A  ) ,  of 
the  water  plus  its  constituents  (Gordon,  et  al.,  1983)  .  The 
contributions  from  the  individual  constituents  can  be  summed 
to  provide  a  total  value  for  each  optical  coefficient. 
Moreover  the  contributions  to  B(A),  b(A),  and  a  (A)  arising 
from  phytoplankton  and  their  pigments  are  assumed  to  be  pro¬ 
portional  tc  chlorophyll  concentration  C.  Taking  a  ratio  of 
L  at  two  different  wavelengths  and  applying  the  as..  ■  ption 

L  (A)  a  3  (A  )b  (A)  /a (A)  ,  we  obtain 
w 
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LwOi)  =  B(.\I)b(:,1)a(x2)  (9) 

L  w  (  X  2  )  B(A2)b(X2_>  a(Xi) 

Because  of  the  non-linearities  involve!  in  the  individual 
constituent  contributions  to  B(X),  b  (A ) ,  and  a  (A  )  ,  we  appeal 
to  a  heuristic  observation  that 

L«Ul>  -  R(C,K,  .  .  .)  (10) 

lu(J2) 

i.e.,  the  ratio  of  two  upwelled  radiances  is  a  function  R 
of  the  chlorophyll  concentration,  C,  and  the  diffuse  attenu¬ 
ation  coefficient,  k,  as  well  as  ether  optically  important 
constituents  of  seawater.  It  was  then  assumed  that  3  is 
related  to  C  through  a  log  linear  model  of  the  form 

Log  C  *  Log  Aq  +  Ai  Log  R(X1,X2)  ,  (11) 

which  was  empirically  fit  to  observations  to  determine  coef¬ 
ficients  Aq  and  Ai .  Thus,  pigment  concentrations  are  com¬ 
puted  from  CZCS  data  using  the  equation 

Ai 

C  -  A  R  (12) 

o  '  ' 

The  empirical  coefficients  presented  by  Gordon, 
at  al.  (1983)  have  been  adopted  by  NASA  and  are: 

Case  I:  C  <  1.5  for  R  (443,550) 


A 


1.  1297959 


Case  II: 


C  >  1.5  for 


F  (443,550)  ,  but 


C  <  1.  5  for  E  (520,  5  50) 


A  =  1.1297950 


A,  =  -1.71 


Case  III:  C  >  1.5  for  E (520,  550) 


A  =  3.3265955 


Ai  =  -2.44 


where  C  is  in  mg/m  . 

b.  Diffuse  Attenuation  Coefficient 

A  similar  development  of  the  algorithm  for  the 
determination  of  the  diffuse  attenuation  coefficient,  k(A), 
is  given  by  Austin  (1981)  .  Like  the  chlorophyll  concentra¬ 
tion  algorithm,  this  algorithm  derives  a  value  based  on  the 

ratio  of  L  at  two  wavelengths.  k  can  be  defined  as 
w 


k(A)  = 


dF ( A , z ) 


F( A  ,  z)  dz 


Equation  (13)  can  be  solved  for  irradianoe  F(a,z)  to  obtain 


F  (  A  ,  z2  )  =  F ( A  ,  z i )  exp  ^ -k( A) (z2-z i ) J 


Hence 


1  ("  F  (  A  ,  z  2  )  i 

k (  A )  =■  - - —  In  , 

2  2“  2  l  J^FCAjZi)  J 
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Empirically  derived  coefficients  from  spectral  data  yieli 


k ( 4  9  0 )  =  0.0833 

(Austin,  1981). 

D.  SION  A I  FACTO  PS 


L  (443) 

w _ 

L  (550) 
w 


.491 

+  0 . 02. 2  m 


(  19) 


Many  factors  have  been  accounted  for  with  these  algor¬ 
ithms  by  either  mathematical  and  empirical  models  or  heuris¬ 
tic  assumptions.  The  determination  of  the  total  radiance 
values  in  the  first  four  channels  of  the  CZCS  allows  us  to 
apply  the  corrections  to  determine  upwelled  radiance.  The 
constituents  of  the  water  which  affect  its  absorption  and 
scattering  properties  are  then  empirically  derived. 

The  distribution  of  phytoplankton  is  controlled  by  many 
local,  mesoscale  and  global  factors,  including  solar  radia¬ 
tion,  global  weather  patterns,  and  ocean  circulation  pat¬ 
terns.  The  mesoscale  events  of  upwellirg  or  eddy 
circulation  can  have  important  regional  effects.  These  fac¬ 
tors  are  too  numerous  and  varied  to  be  modelled  on  a  theo¬ 
retical  basis.  However,  empirical  modelling  can  produce 
relatively  accurate  and  consistent  results. 

The  measurement  of  these  bio-optical  parameters  from 
space  allows  us  to  remotely  determine  their  relationships  to 
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physical  everts  in  the  regions  under  stuiy.  Tine  seals-, 
spatial  scales  and  specific  features  can  be  discerned  us  in a 
the  known  (or  hypothesized)  relationships  between  inherent 
optical  properties  of  the  ocean  water  constituents  and  the 
forcing  involved  in  their  distribution. 


f  S 


IV 


EXFIrvICAl  OS  T  r.O  A  C \  L  ^  u.iCl*  j.i 

- TTT770IT5 - 


A.  INTHO  DUCT  ION 

The  concept  of  principal  component  analysis  has  been 
presented  and  utilized  in  different  forms  over  the  past 
eighty  years.  Fitting  a  line  to  a  data  set  was  usually 
accomplished  using  a  least  squares  method.  Distances  to 
this  line  from  each  point  were  measured  parallel  to  an  arbi¬ 
trarily  set  axis.  From  the  early  work  of  Pearson  ( 190 1)  , 
this  method  was  adapted  sc  that  the  perpendicular  distances 
from  each  point  to  the  best  fit  line  were  measured.  Figure 
4  illustrates  this  difference  and  shows  that  the  first 
method  is  tied  to  a  cocrdinate  system  while  the  Pearson 
approach  is  independent  of  coordinate  systems.  This  new 
method  laid  the  foundation  for  the  development  of  principal 
component  decomposition  techniques.  These  techniques  have 
since  been  utilized  in  many  forms  and  referred  to  by  similar 
names  in  a  number  of  disciplines.  Applications  in  psychol¬ 
ogy  by  Eckert  and  Young  (  1936 , 1  93 9)  ,  although  somewhat  dif¬ 
ferent  in  their  development,  contain  the  essential  elements 
of  data  analysis  and  principal  component  decomposition  as 
used  in  geophysical  disciplines  today. 


Meteorological  applications  by  Lorenz  (1956) ,  Kurz- 
bach  (1967)  ,  and  Rinne,  et  al.  (19  79)  demonstrated  the  ccnvia- 
nance  of  representing  a  large  climatological  field  with  a 
smaller  set  of  values.  These  areas  include: 

(1)  Non-linear  statistical  prediction  (Lorenz,  1956), 


(2)  Ncn-linear  dynamical  prediction  (Lorenz ,  1956)  , 


(3)  500  mb  height  field  representation  (Rinns,  et  al., 

1979) ,  and 


(4)  Sea  level  pressure,  surface  temperature,  and  pre¬ 
cipitation  pattern  representations  (Kutzbach,  1967) 
Other  uses  of  EOF  analysis  techniques  in  oceanography 
include  the  representation  of  ocean  color  spectra  (Mueller, 
1976)  and  of  wave  spectra  (Aranu vachupun  and  Thorton,  1983)  . 

The  principal  difficulties  encountered  in  principal  com¬ 
ponent  analysis  problems  relates  to  the  selection  of  the 
'meaningful*  subset  of  components  and  to  their  physical 
interpretation.  Methods  of  selection  of  the  principal  com¬ 
ponents  are  alsc  widely  varied.  Preisendorf er ,  et  al.  (1981) 
discussed  two  methods  which  together  involve  seventeen  dif¬ 
ferent  testing  rules.  Empirical  selection  of  a  cutoff  value 


for  variance  or  forcing  factors  can  also  be  utilized. 


Visual  ir. spection 


of  the  data  which  laads  to  a  clear  cut 


(albeit  subjective)  choice  is  also  an  option. 

The  enormous  data  volume  inherent  in  satellite  data  sets 
begs  application  of  the  techniques  of  principal  component 
analysis.  Principal  component  analysis  techniques  often 
allow  the  efficient  representation  of  a  large  data  set  by 
its  first  few  principal  components  with  a  negligible  loss  of 
information.  The  advantage  gained  is  reduction  in  the  num¬ 
ber  of  variables  needed  to  represent  the  data.  Reducing  a 
data  set  to  its  principal  components  can  also  aid  in  the 
interpretation  of  the  data  by  separating  noise  from  the  sig¬ 
nal.  Principal  component  analysis  theory  can  be  applied  to 
preliminary  explorations  within  a  relatively  unstructured 
domain  of  knowledge,  one  in  which  the  fundamental  laws  gov¬ 
erning  the  processes  under  study  are  still  being 
defined.  (Preisendor f er,  et  al.,  1981) 

A  brief  review  of  the  EOF  analysis  follows  to  provide 
background  for  the  later  analyses.  The  reader  is  referred 
to  Priesendorfer,  et  al.  (1981)  for  a  more  complete  devel¬ 
opment  and  history.  The  following  matrix  algebra  notation 
is  adopted  throughout  this  thesis. 

1.  Nc  underscore  denotes  a  scalar . X 

2.  A  straight  line  underscore  denotes  a  vector.... £ 
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3. 

4. 


A 

A 


curved  line  unde 
straight  line  cv 


core  denotes  a  matrix . . 

oar  denotes  a  ;nean  value... ? 


5.  The  use  of  a  superscript  "T"  denotes  a  matrix 
transpose. 


B.  EOF  EQUATIONS 

1.  Haw  Data  Conversion 

Following  Preisen  dorf  er ,  et  al.  (1981)  let  F'  be  the 
raw  ( unce ntered)  data  matrix. 


F' 


f ' (1,1)  f ' (1,2)  ...  f ' (l,p) 
f ' (2, 1)  f ' (2,2)  ...  f ' (2,p) 


(n,l) 


f ’ (n,2) 


f ' (n , p) 


(17) 


where  f'  (i,  j)  is  the  measurement  in  the  i'th  time  point  and 
j'th  spatial  position.  In  the  present  investigation,  each 
member  of  will  correspond  to  an  optical  parameter  meas¬ 

ured  ty  the  CZCS  at  a  particular  time  and  spatial  position. 
To  convert  the  raw  data  matrix,  F',  to  a  centered  data 
matrix,  F,  the  temporal  means  are  computed  and  subtracted 
frcm^F'.  The  temporal  mean  vector  f(x)  is  calculated  as 

—  V'  (18) 

f  (x)  =  n  =  l  f ' (t  ,x) 


The  centered  data  matrix,  F.  is  then  defined  as 


f '  (1,1)  -  f ( 1 )  f ' (1,2)  -  f(2)...f'(l,p)  -  f  ( ?  }  | 

f  '  ( 2  ,  i )  -  ill)  f  ’  (2,2)  -  £(2)  (2, p)  -  ! 


L  f  '  (n  ,  1 )  -  f(l)  f  '  (n  ,  2)  -  f  (2)  .  .  .  f  '  (n,p)  -  f(p)J  • 

Each  element  f(t,x)  of  J  consists  of  a  raw  data  measurement 
with  the  temporal  mean  removed.  The  centered  data  matrix  , 
F,  can  be  written  as 


F  =  F' 


(20) 


where  ?  is  the  matrix  containing  as  rows  the  transpose  of 
the  mean  vector  fT=  £  f  (1 )  ,  f  (2) , .  . . ,  f  (p)  j  . 

2.  Principal  Direction  of  Scatter 

Tc  find  the  direction,  £i,  (in  the  space  domain) 
along  which  the  scatter  (or  variance)  of  the  data  set  is 
greatest,  consider  the  projection  of  the  data  vectors  fjt) 
along  an  arbitrary  direction  ej 

D  (  t  ,  e  ! )  =  f  T(  t )  e  J 


(21) 


- ...  •  • 


Sg ua  ring  this  length  and  summing  over  all  r.  observan icr.s. 


41  Via  a  ul  te  c  S  -i  £  (i  C  £  oCdtTT'si 


.  i  -  i  i  iO  1‘.  .j 


tior. ,  e  i #  namely 


D2(e,)  -  £  i  [**<«>.,  j 


(22a) 


The  righthana  side  of  (22a)  can  be  expanded  to  yield 
n 


D  ( e  j )  = 


=  £  ("  eTf(t)fT(t)e 

t  =  i  L  1  x. 


T 

■  Si 


n  — 

£  f(t)fi(t)  ex 
_  t  =  l" 


(22b) 

(22c) 


The  next  step  is  to  define  the  "Scatter  Matrix",  S, 


S  =  F  F 

r*  /sj  r*> 


(23) 


with  elements 


(i,j)  =  £  f(OfT(t) 
t  =  i 


(24) 


Expanding  the  above  equation  produces  for  each  member  of  J5 
n 

T 


(i,j)  =  £  (f • (t,i)-f (i)) (f ' (j))  .  (25) 

t  =  1 


If  the  matrix  is  normalized  by  dividing  by  (p  -  1),  then 

when  i/  1,  the  members  of  S  are  covariance  values,  and  when 
i  =  j  (the  trace  of  the  matrix)  the  members  are  variance 
values  (i.e.  each  element  is  the  variance  of  f  at  a  single 


The  scatter  matrix,  S, 

la 
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spatial  grid  point)  . 


is  symmetric. 


Therefore,  it  generally  has  p  no  r.-n  e  gative  eigenvalue-  i .  .  'g 
=  1,  p)  and  associated  eigenvectors  Sj  (j  -  •  ••»  p)  / 

o-cvided  that  the  rank  of  S  is  equal  to  ?. 

3.  principal  Com£ onent,  Eig envalue  and  Eigenvector 
82££§sen^afcn  ~  "* 

The  first  principal  component  or  an  observation  vec¬ 
tor  f  is  defined  to  be  the  linear  combination 


a  1  -  euf  i  +  e  1 2  f  2  +  •  •  •+  elpfp  =  ®  1 


(26) 


whose  sample  variance 
p  P 


(27) 


Sa  =  ^2  S  6 1  ie  1  j  Si j  =  -U-1  . 

1  i  -  1  j  -  1 

is  a  maximum  for  all  possible  vectors  e,  subject  to  the  con¬ 
straint  that 


®  1  ®  l 


=  1 


(28) 


Introducing  the  Lagrange  multiplier  Xj,  the  maximum  variance 
must  satisfy 


3 

5e7 


+ 

e  1 


X 1 ( 1-efei 
=  0 


) 


(29) 


For  ncn-trivial  solutions 


\x,  must  be  chosen  such  that 


S  -  xa 


o 


(30) 


and,  therefore.  ,  is  an  eigenvalue  of  S,  and  Sj  is  ids 
associated  eiaer. vector.  Furthermore 


s 

a, 


Ei 


■l  Ei 


and  since  e  ie  i  =  1 


T 

e  S  e  =  A 
"1  J  ~1  1 


(31) 


(32) 


i .  e . ,  the  first  eigenvalue  of  is  interpretable  as  the  sam¬ 
ple  variance  of  S.  If  we  expand  this  development  to  the 
other  eigenvalues  and  eigenvectors  of  S,  we  obtain 


S  =  [-j]  ;  j  “  1’*‘*P  ’ 

the  eigenvector  matrix  and 


(33) 


L  -  D[X.];  J  -  !.-»•  ’  <3»> 

the  diagonal  eigenvalue  matrix. 

Thus , 

i  i  -  s  t  •  <35> 

T  T 

In  terms  of  E  the  constraint  e,  e  =  1  becomes  E  E  *  I.  Where 

</>  —i—i  >/'>/'  <✓> 

I  is  the  identity  matrix.  Therefore,  if  we  multiply  both 

O' 

T 

sides  of  equation  (35)  by  £  we  obtain 

l  ’  I  k  sT  • 


(36) 


Now,  using  the  definition  provided  for  f  and  equation  (20)  , 


the  principal  component  matrix  can  be  defined  as 

ll  .  (37) 

This  is  the  desired  principal  decomposition  of  F  where 

^  k  =  l  l  <38a> 

£  £  £  £  (38b) 

*  $T  $  E  (38c) 

*  f  l  \  (38d) 

*  k  *  (38e) 


C.  PARTITIONED  EOF  ANALYSIS 
1 .  Purpose 

The  EOF  analysis  method  outlined  above  works  very 
well  for  a  large  continuous  data  set.  However,  geophysical 
data  sets  are  rarely  continuous.  In  the  case  of  satellite 
data,  cloud  cover  results  in  many  gaps.  Sometimes  these 
gaps  can  he  bridged  by  linear  interpolation,  e.g.,  when  they 
are  small  and  surrounded  by  good  data.  Often  this  is  not 
the  case  and  so  a  scheme  of  utilizing  non-contin uous  data  is 
necessary.  Here,  the  purposes  of  partitioning  are: 


57 


(1)  To  maximize  the  sample  size  in  the  presence  of 
cloud  cover,  thus  allowing  statistical  computations  for 
subregions ; 

(2)  To  highlight  spatial  structures  of  variance  fea¬ 
tures  locally,  before  absorbing  them  into  the  modes  of 
the  overall  domain;  and 

(3)  To  achieve  computational  convience. 

Briefly,  partitioning  permits  EOF  analysis  using 
small  subsets  of  the  overall  data  set.  These  subsets  are 
partitioned  to  yield  continuous  data  in  each  subdomain.  An 
EOF  analysis  is  completed  on  each  individual  subset,  and  an 
eigenvalue  matrix,  an  eigenvector  matrix  and  a  principal 
component  matrix  are  obtained.  The  next  step  is  to  perform 
an  EOF  analysis  to  join  the  principal  components  of  the  sub¬ 
sets.  This  second  EOF  analysis  produces  'joining  functions' 
which  relate  two  non-overlapping  subsets. 

2.  Rules  and  a e t ho d s 

When  performing  the  partition  of  any  data  set  cer¬ 
tain  rules  must  be  observed  tc  maintain  the  statistical 
reliability  cf  the  computations.  Two  obvious  and  basic 


rules  are; 


(1)  The  minimum  partition  siz a  (number  of  pixels)  ?.us* 
cs  greater  than  or  equal  to  the  sample  size  (i.e.,  if 
there  are  25  sample  days  each  partition  must  have  25  or 
more  pixels).  In  practice,  the  spatial  dimension  will 
fce  required  to  be  significantly  greater  than  the  sample 
size . 

(2)  The  partition  size  should  net  be  so  small  that  the 
spatial  structure  is  dominated  by  noise  (e.g.  ,  a  parti¬ 
tion  boundary  will  not  be  placed  in  a  major  feature, 
such  as  a  front  or  eddy  of  length  scale  much  less  than 
the  partition  size). 

The  methods  involved  in  the  partitioning  are  subject 
to  the  abeve  principles,  together  with  a  general  understand¬ 
ing  cf  the  physical  processes  occurring  in  the  study  domain. 
Four  tracklines  at  35N,  35  22’N,  35  40’N  and  35  53'N,  were 
used  to  aid  in  this  initial  trial  of  partitioning  (Fig.  5)  . 
The  radiance  values  and  computed  optical  parameters  along 
each  trackline  were  plotted  versus  distance  from  the  coast. 


These  plots  were  aligned  to  pictcrially  represent  the  data 
and  its  gaps  (due  to  clouds).  The  partitioning  scheme  was 


then  applied  tc  try  and  produce  subsets  that  were  as  ccji- 
pieze  (continuous)  as  possible  over  the  time  domain.  The 
total  data  matrix  is  thus  partitioned  inzo  ?  subdomains 


It  "  £p  ;  p  =  1  ’  '  * ‘  ,P  ’ 


where  the  subscript  t  denotes  the  zozai  data  set  and  sub¬ 
script  p  denotes  the  partitions  of  the  data  set.  Bach  F  is 
zhe  data  matrix  for  grid  points  falling  within  grid  parti¬ 
tion  p,  and  contains  all  time  points  for  which  complete  data 
ware  acquired  in  that  subdomain. 

3.  Equation  Development 

The  EOF  decomposition  discussed  in  B,  is  applied  to 
each  parxiticn  separately,  such  that  for  each  partition,  p, 
zhe  scatzer  matrix  is  given  by 


lp 


and  from  (36) 


S  =  E  L  E  (41) 

'Vp  a>p  p  a.p  f 

where  jS  *  Eigenvector  matrix  of  the  spatial  partition  p. 


p  *  Eigenvalue  matrix  for  partition  p. 


The  matrix  of  principal  components  for  each  partition  is 
given  by 


A  -  F  E 

O/p  'Xi  p 

Equations  (38)  require  that 


h 


T 

-  A  A 

p  'X,  p  %  p 


and  so  (4  1)  can  be  written  as 


T  T 

S  =  E  A  A  E 

■tp  ^jp  ^  p  ^p  'Vp 


(42) 


(«3) 


(44) 


Relating  this  to  a  global  scatter  matrix,  S-, 

\S)  li 


(SYM) 


c7  s 

'\j*-  p  a»p 


(45) 


Where  1 2  represents  the  matrix  of  covariances  between  grid 
points  in  domains  1  and  2,  and  so  forth. 

Now  the  joining  process  is  developed.  For  any  num¬ 
ber  cf  partitions  (two  are  used  in  this  development)  , 


Si  -  |!  *  £1  $! 

and 

$2  3  £2  £2  *  $2  $2  $2  ^2 


(46) 


(47) 


s 


s 


where  the  subscripts  denote  partition  number. 

For  this  combined  set  the  scatter  matrix  is  giver,  by 


S 


1  0 


T  T 

F i  F i  F2  Fi 

fl  %z  fz  lz 


(48) 


Using  equations  (46)  and  (47)  , 


’$1  #  ‘ 

>1  $1  £  fa 

T  T 

$  $2 

$1  £2  £2  £2 

(49) 


The  joining  functions,  are  defined  as  the  eigenvectors  of 
the  central  matrix  given  in  (49)  .  Finally  using  (4  1)  ,  (43) 
and  (44)  ,  (49)  becomes 
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The  joining  functions,  ^J,  relate  the  separate  subscenes  to 
each  ether  across  an  overall  study  domain.  The  interpreta¬ 
tion  of  these  functions  should  allow  examination  of  the  var¬ 
iations  that  occur  throughout  the  domain  and  localized 
effects  on  the  individual  partitions  of  the  domain. 
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D.  INTERPRETATION 

The  coordinate  system  defined  by  the  eigenvectors  gives 
the  domain  for  the  principal  components.  In  the  present 
study,  as  well  as  most  geophysical  applications,  the  princi¬ 
pal  components  can  be  thought  of  as  temporal  amplitudes  and 
the  eigenvectors  as  their  special  modulators  (Preisendorf er, 
ei  al.,  1981) . 

The  i-th  principal  component  is  that  linear  combination 
of  the  data  field  which  explains  the  i-th  largest  portion  of 
the  total  field  variance.  Essentially  the  eigenvectors 
define  a  direction  of  variance,  while  the  principal  compo¬ 
nents  give  the  amplitude  of  the  variance  in  the  direction  of 
the  associated  eigenvector. 

Once  a  data  set  has  teen  reduced  to  a  set  of  eigenvec¬ 
tors  and  associated  amplitudes,  the  question  of  signal  ver¬ 
sus  noise  arises.  A  decision  as  to  which  components  of  the 
data  field  have  significance,  and  which  components  of  the 
data  field  have  no  physical  meaning  must  be  made.  Seme  sort 
of  a  selection  process  must  be  defined  and  applied.  Bases 
for  these  selection  processes  should  have  their  roots  in  the 
physical  processes  being  studied. 


V.  RESULTS 

A.  INTRODUCTION 

The  major  focus  of  this  thesis  was  to  achieve  a  first 
step  towards  the  analysis  cf  the  obtained  data  set.  Much  of 
the  energy  in  producing  these  results  was  directed  toward 
the  processing  of  the  data  to  a  usable  form  for  EOF  analy¬ 
sis.  Appendices  A  and  B  give  a  detailed  discussion  of  the 
processing  techniques  utilized  and  an  accounting  of  all 
adjustments  applied  to  the  data. 

The  results  obtained  in  this  thesis  encompass  three  dis¬ 
tinct  areas.  The  first  result  emerged  from  the  data 
processing  and  a  discovery  of  the  breakdown  in  the  black 
ocean  assumption.  The  remaining  areas  are  interrelated  as 
one  deals  with  the  data  set  prior  to  the  EOF  analysis,  while 
the  ether  attempts  to  relate  this  to  a  statistical  meaning 
using  EOF  analysis  methods. 

B.  CORRECTIONS  FOR  NON-ZERO  L  (670)  IN  COASTAL  WATERS 

w 

Preliminary  examination  of  this  data  set  showed  that, 
near  the  California  coast,  the  assumption  that  Lw(670)  *  0 

breaks  down  (Chapter  III,  Section  C)  .  This  finding  pre¬ 
sented  a  need  for  an  adjustment  algorithm. 
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k t  pixels  where  the  upwelled  radiance  at  670  r.s  , 
L  (670),  is  significantly  greater  than  zero,  calculated  val- 

W 

2 

ues  cf  L  (443)  are  often  <  0.01  (mw/(c:n  -sr-m) )  (apprcxi- 

W 

mafsly  1  digital  ccunt  in  CZCS  channel  1).  This  is 
unreasonable  even  in  moderately  turbid  ocean  waters.  Smith 
and  Wilson  (1981)  observed  that  in  coastal  waters  off  Cali¬ 
fornia,  where  pigments  and/or  sediment  concentration  are 
relatively  high,  it  is  not  uncommon  for  the  subsurface 
upwelling  radiance  L w (6 7 0 )  to  be  non-zero.  They  developed 
an  iterative  procedure  to  account  for  this,  which  is  similar 
to  that  developed  independently  and  used  in  the  present 
processing . 

The  procedure  involves  two  major  steps. 

In  the  first  step,  which  is  invoked  when  L  (443)  <  0.01: 

W 

1.  Set  1^(443)  =  0.01,  a  minimal  value  for  daylight 

tackscatter  and  slightly  less  than  one  digital  count  in 
CZCS  channel  1. 

2.  Decrease  L  a(670)  and  increase  Lw(670)  to  be  consis¬ 
tent  with  the  new  value  of  Lw(443)  (using  equations  (4) 
and  (1)). 

3.  Recalculate  L  ^( 5 2 C)  and  Lw(550). 

4.  Recalculate  C  u  C2,  K(4  90),  K(520). 
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he  second  step  is  based  on  rhe  assumption  that  the 


C  2  algorithm  (eguation  (12),  case  III)  is  robust  and  insen¬ 
sitive  to  moderate  errors  in  L  ^ (67 0)  .  This  assumption  was 
supported  by  sensitivity  calculations  which  showed  Cz  values 
to  vary  by  less  than  30  percent  for  wide  variations  in  e  (X 
,670).  If  Ci  -  C2  >  0.5  and  C  i  >  2,  tha  correction  is  pre¬ 
sumed  unreasonable  due  to  L  (670)  being  too  large  (and  cor- 

a 

respcndingly ,  L  (670)  being  too  small).  Next,  values  of 

W 

K (490)  and  K  (520)  are  estimated,  which  are  consistent 
wi  th  C  2: 

1.  Estimate  a  ratio  (4  43) /L ^(550)  consistent  with 

C  2  by  inverting  the  C  i  algorithm  (Equation  ( 1 2)  with  case 
I  coefficients)  using  the  C  2  values. 

2.  Increase  L  (670)  and  decrease  L  (670)  to  be  consis- 

W  B 

tent  with  the  new  values  of  Lw(443). 

3.  Recalculate  L  (443),  L  (520),  and  L  (550) 

w  w  w 

4.  Recalculate  C  1#  C2,  K(4  90),  K(520). 

5.  Iterate  this  procedure  until  Cj  and  C2  agree. 

Data  acquired  aboard  the  R/V  Acania  during  the  Optical 
Dynamics  Experiment  (ODEX)  provide  a  tentative  basis  for 
assessing  the  validity  and  performance  of  the  above  adjust¬ 
ment  algorithm.  In  Figure  6  values  of  1/K(490)  calculated 
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from  C2CS  data,  acquired  cn  16  October  1  982,  are  compares 


with  preliminary  calculations  of  1/K(490)  from  selected  CD2X 


stations.  The  transect  shewn  is  along  35N  (partition  4). 


Stations  24  and  25  were  occupied  1.5  hours  before  the  Nim¬ 


bus-7  CZCS  observation,  and  2  hours  after  it,  respectively. 


Station  21  was  occupied  9  hours, and  stations  19  and  13  one 


and  two  days,  respectively,  prior  to  the  satellite  pass. 


The  1/K(490)  values  at  these  stations  were  calculated  from 


the  graphical  displays  of  raw  irradiance  profiles  (at  a 


wavelength  of  490  nm)  presented  in  the  preliminary  R/V 


ACANIA  0 0 EX  CRUISE  REPORT  (Mueller,  Zaueveld  and  Smith 


1982)  . 


Panel  6a  compares  the  CZCS  and  in  situ  1/K  (490)  values 


before  the  above  adjustment  was  applied,  and  figure  6b  com¬ 


pares  them  after  the  correction.  Agreement  in  both  cases  is 


excellent  in  the  transparent  waters  at  stations  21,  24  and 


25:  no  adjustment  for  1^(670)  was  required  in  this  region. 


In  the  inshore  portion  of  the  transect,  however,  agreement 


is  obviously  poor  before  the  L  (670)  adjustment,  and  much 

W 


improved  afterwards.  This  result  is  preliminary,  and  sub¬ 


ject  to  possible  revision  by  cognizant  ODEX  investigators 


when  their  data  have  been  brought  to  publishable  form. 
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Na verthaless,  the  L  (570)  adjustment  algorithm  sc  overwhelm- 

W 

ingly  improves  the  CZCS  estimates  of  K(490)  that  its  use  in 
this  thesis  project  is  fully  justified  and  essential. 


500  400  300  200  100  0 


DISTANCE  (KM) 


.  Comparison  Plots  For  1/K  (490)  Between  Track  4  and 
Selected  ODEX  Stations.  Panels  (a)  and  (b) 
Respectively  are  Before  and  After  Adjustments  for 
Non-Zero  Values  of  L  (670)  . 
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Figure  6 
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C.  DATA  STRUCTURE 

Figures  7  through  10  shew  the  optical  depth  parameter, 
1/X{490)  =  Z9  0  (490)  along  each  track  for  the  available  data 
scenes.  (Gordon  and  acCluney  (  1975)  showed  that  Z90(x)  is 
the  depth  over  which  90  percent  cf  L w ( A)  is  backscattered. ) 
The  plots  are  oriented  sc  that  the  coast  is  on  the  right- 
hand  side  (positive  x)  ,  while  time  of  the  data  scene  gees 
from  earliest  to  latest  in  the  positive  direction  along  the 
ordinate  cf  each  figure.  The  scale  of  1/K  (490)  is  in 
meters. 

Chapter  II  and  III  give  background  into  the  oceanography 
of  the  region  and  how  that  can  be  related  to  ocean  optical 
parameters.  The  structures  depicted  in  figures  7  through  10 
will  be  discussed  in  terms  cf  ocean  eddy  and  front  visuali¬ 
zations  which  result  from  these  relationships.  Relatively 
high  values  of  1/K(490)  indicate  water  with  lower  concentra¬ 
tions  of  chlorophyll  and  sediment.  In  general  these  concen¬ 
trations  may  be  expected  to  decrease  with  distance  offshore. 
Abrupt  changes  in  1/K  (490  )  are  usually  associated  with  ocean 
frontal  structure  and  eddies. 


0 
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P  a-t  iljLon  \  (Zonal  Transect  at  35  53 N) 

The  only  data  available  from  1979  (23  Nov) 


troa 


winter  and  shows  relatively  little  structure,  (Fig*  7)  • 
This  image  was  obtained  after  the  end  of  the  upwelling  sea¬ 
son  and  the  surface  waters  were  homogeneous  to  at  least  225 
km  offshore. 

The  1980  data  series  shows  more  structure.  Begin¬ 
ning  on  17  May  80,  an  eddy  of  approximately  40  km  diameter 
was  centered  approximately  180  km  offshore.  Sixteen  days 
later  the  entire  track  shows  several  eddy-like  features 
ranging  in  size  from  4  km  to  20  km.  Three  days  later,  on  6 
June  80,  the  track  has  lost  much  of  this  structure,  although 
1/K  (490)  generally  increases  in  the  offshore  direction. 
This  trend  persisted  and  strengthened  slightly  through  June 
1980.  By  1  August  1980,  a  distinctive  pattern  had  developed 
with  nearly  uniform  turbid  waters  adjacent  to  the  coast,  and 
an  abrupt  (15  km)  frontal  transition  to  much  more  transpar¬ 
ent  waters  at  a  distance  approximately  95  km  offshore.  This 


pattern  is  suggestive  of  the  zonal  seals  of  bio-optical 
response  to  coastal  upwelling  over  a  single  season. 
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The  1982  scenes  are  ail  in  the  lane  fall  to  early 
winter.  The  transect  from  30  September  1982  shows  uniformly 
turbid  water  (5  to  10  m  optical  depth)  with  little  struc¬ 
ture,  except  for  a  weak,  dear-water  eddy  signature  200  km 
offshore.  Between  30  September  and  16  October,  more  trans¬ 
parent  (20  to  25  m  optical  depth)  water  intruded  from  off¬ 
shore  to  within  190  km  of  the  coast,  with  a  frontal 
transition  region  of  scale  approximately  15  km.  Over  the 
ensuing  month,  this  clear  intrusion  appeared  to  evolve  into 
a  field  of  less  organized,  eddy-like  anomalies  with  scales 
ranging  frcm  10  to  50  km.  The  2-dimensional  character  of 
this  pattern  evolution  should  be  studied  in  a  future  analy¬ 
sis  cf  this  data  set.  Such  an  analysis,  which  is  beyond  the 
scope  cf  the  present  thesis,  may  indicate  whether  these 
changes  are  test  interpreted  as  breakdown  of  a  spatially 
continuous  intrusion  of  offshore  water,  or  as  simply  due  to 
advection  transporting  entirely  different  water  mass 
features  into  view  at  this  location. 

2.  Partition  2  (Zonal  Transect  at  35  40  N) 

Track  two  (Fig.  8)  is  located  along  latitude  35 
40’  N  (south  of  the  first  track)  . 


Again,  the  1979  data  set  is  limited  to  one  Say  ia*:  = 
in  the  year  and  is  virtually  featureless. 

The  1980  series  starts  in  late  spring  (6  May  1980) 
and  progresses  through  late  summer  (1  August  1980)  .  The 
same  general  features  found  in  the  northern  track  are  also 
evident  in  track  2.  However,  the  general  appearance  of  the 
frontal  boundary  is  broken  by  additional  eddy  structures. 
In  addition,  organized  ocean  frontal  structures,  which  are 
apparently  related  to  similar  features  in  Partition  1  data, 
are  here  displaced  approximately  15  km  further  offshore. 

The  1982  scenes  from  partition  2  show  the  same 
front/eddy  development  presented  in  the  data  from  the  north¬ 
ern  transect  (partition  1).  An  offshore  frontal  boundary 
formed  between  30  September  1982  and  16  October  1982,  and 
evolved  to  a  less  organized  pattern  of  eddy  signatures. 

3-  Partition  3  (Zonal  Transect  at  35  22 N) 

Along  track  3  (Fig.  9)  the  features  already  dis¬ 
cussed  for  the  previous  partitions  occurred.  An  additional 
feature  here  is  a  strong  clear-water  eddy  signature  approxi¬ 
mately  220  km  offshore,  evident  on  the  third  of  June  1980. 
The  data  suggest  an  eddy-like  intrusion  of  transparent  off¬ 
shore  water,  approximately  40  km  in  diameter.  (A  feature  of 


similar  diameter,  but  centered  approximately  40  km  clorer 
inshore  appears  simultaneously  in  partition  4  (see  b--low))  - 
This  feature  is  short-lived  ,  however,  as  no  evidence  of  it 
remains  three  days  later.  This  strange  phenomenon  is 
another  candidate  for  future  investigation  in  a 
2-dimensional  analysis. 

In  general,  frontal  boundaries  along  partition  three 
are  farther  offshore  in  both  the  1980  and  1982  series  than 
apparently  related  frontal  signatures  in  partitions  one  and 
two.  Again,  this  may  suggest  the  importance  of  bottom 
steering  of  the  mean  flow  over  diverging  isobaths. 

4.  Partition  4  (Zonal  Transect  at  35  00 N) 

The  southernmost  track  at  35  00  N  (Fig.  10)  has  sim¬ 
ilar  patterns,  including  the  offshore  displacement  of  fron¬ 
tal  boundaries.  The  location  of  the  frontal  boundary  is  at 
its  farthest  offshore  position  in  this  track. 

The  transient  eddy  feature  of  3  June  1980,  as  dis¬ 
cussed  in  partition  3,  appears  here  to  have  a  latitudinal 
extent  large  enough  to  span  partitions  3  and  4  which  are 
roughly  45  km  apart.  This  is  consistent  with  the  zonal 
scale  of  the  feature  (approximately  40  km). 
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There  also  is  a  longitudinal  displacement  of  approximately 
40  km  which  suggests  an  eddy  of  slightly  oblong  shape  that 
roughly  parallels  ths  coast.  Questions  of  whether  such  fea¬ 
tures  are  coherent  and  continuous  from  one  transect  to 
another  can  fce  easily  resolved  in  future  2-diaen  sior.al 
analyses. 


It  is  reasonable  to  expect  a  band  of  surface  water 
with  low  optical  depth  to  lie  adjacent  to  the  coast,  due  to 
enhanced  nearshore  biological  productivity  during  the 
upwelling  season  and  due  tc  northward  transport  of  more  tur¬ 
bid  water  masses  by  the  Cavidson  Inshore  Current  in  winter. 
In  contrast,  offshore  waters  tend  to  be  far  more  transpar¬ 
ent,  at  least  during  the  seasons  covered  by  this  data  set. 
It  may  be  anticipated  therefore,  that  an  obvious  optical 
front  will  persistently  delineate  the  boundary  between  what 
may  be  classified  as  nearshore  and  offshore  bio-optical 
regimes,  and  that  intrusions  of  eddy- like  surface  water  fea¬ 
tures  from  one  regime  to  the  other  will  be  illuminated  by 
optical  contrast.  Additionally,  horizontal  gradient  struc¬ 
ture  in  bio-optical  processes  will,  within  each  regime, 
often  (but  not  always)  accompany  the  physical  structure 
associated  with  <cean  fronts  and  eddies  and  produce  optical 


depth  structure  cf  similar  scales.  The  scales  cf  previously 
observed  structures  in  tie  California  Current  region  may 
thus  be  expected  to  be  present  in  the  horizontal  structure 
of  optical  depth.  It  is  clear  from  the  above  discussion  and 
the  review  cf  Chapter  II  that  this  is  in  fact  the  case.  In 
the  next  section  these  structures  will  be  discussed  frcm  a 
statistical  viewpoint  using  EOF  analysis. 

D.  EOF  ANALYSIS 

1  •  Eigenvalues  and  Degrees  of  Freedom 

Table  II  is  a  listing  cf  the  eigenvalues  for  the 
spatial  covariance  matrix  calculated  from  the  data  in  each 
partition  analyzed.  The  cumulative  percentages  cf  the  total 
variance  are  included.  The  first  ten  eigenvalues  are  listed 
here.  In  each  case,  they  account  for  roughly  S8  percent  of 
the  variance.  The  eigenvalues  are  presented  graphically  in 
figures  11  through  1 4,  which  include  both  the  eigenvalues 
listed  in  Table  II,  and  the  additional  eigenvalues  that  con¬ 
tain  the  nciser-appearing,  higher  spatial  frequencies  of 
variability  (and  together  account  for  lass  than  2  percent  of 
the  total  sample  variance).  We  have  assumed  that  98  percent 
of  the  variance  is  an  adequate  cutoff  for  calculations. 
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Ihe  major  features  seen  ir.  the  eigenvalues  are  rhu 
differences  in  the  first  value  between  partitions  1,  2,  and 
3,  and  partition  4.  The  variance  in  partition  4  is  mere 
evenly  distributed  over  the  first  three  eigenvalues.  All 
four  partitions  have  reached  roughly  90  percent  of  the  vari¬ 
ance  by  the  sixth  eignevalue.  The  difference  in  the  struc¬ 
ture  of  the  eigenvalues  suggests  that  partition  four  is 
either  affected  by  additional  factors  not  found  in  the 
northern  partitions,  or  that  some  factors  which  influence 
the  northern  partitions  are  absent  here. 


2.  Data  Reconstruction  Psin g  Eigenvectors  and  Principal 
Components!”  ”  ~  ~ 

Before  the  eigenvectors  and  principal  components  are 
interpreted,  how  they  are  combined  to  reconstruct  a  particu¬ 
lar  observation  is  explained.  Recall  that  each  eigenvector 
defines  a  direction  of  spatial  varibility,  and  that  its 
associated  principal  components  represent  the  amplitude  of 
variations  in  that  direction  at  certain  time  points.  In  the 
present  context  a  "direction"  takes  the  form  of  1/K(490) 
variations  that  are  coupled  at  all  grid  points  of  the 


domain,  and  "direction"  in  this  sense  may  be  best 
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TABLE  II 

Eigenvalue  Data  for  Partitions  1  through  4 


Partit 

ion  1 

Parti 

tion  2 

Eigenvalue 

2 

Cumulative 

Eigenvalue 

2 

Cumulative 

(m  )- 

Percentage 

(m  ) 

Percentage 

271. 40 

4  1.  60 

365.90 

43.73 

174. 90 

68.41 

217.40 

69.67 

72.  17 

79.47 

99.03 

81.49 

37. 72 

85. 25 

4  5. 60 

86.94 

26.  SO 

89.37 

28.85 

90.39 

17.69 

92.08 

20.39 

92.83 

14.80 

94.  35 

15.12 

94.64 

10.59 

95.97 

13.63 

96.27 

7.  02 

97.  05 

7.90 

97.21 

5.08 

97.83 

5.98 

97  .92 

Partition  3 

Parti 

tion  4 

Eigenvalue 

Cumulative 

Eigenval  ue 

2 

Cumulative 

2 

(m  )  ♦ 

P  ercentage 

(m  ) 

Percentage 

557.40 

54.  18 

273.70 

35.52 

210.60 

74.65 

183.60 

59.35 

97.76 

84.  15 

119.80 

74.90 

58.08 

89. 79 

44. 14 

80.63  . 

32.93 

92.99 

40.87 

65.93 

20.72 

95.00 

35.96 

90.60 

13.86 

96.  35 

24.42 

93.77 

9.48 

97. 27 

13.38 

95.51 

8.  29 

98.  08 

8.22 

96.58 

5.74 

98.64 

7.00 

97.49 

illustrated  either  as  a 


curve  (for 


one-dimensional 


transects)  or  contour  plots  (for  two-dimensional  domains). 


Given  a  temporal  mean  value  for  each  grid  point  of  a  domain. 


the  eigenvector  multiplied  by  the  principal  component  for  a 
specific  time  yields  a  modifier  tc  the  mean  signal. 

Por  example,  to  view  the  contributions  of  the  first 


five  eigenvectors  to  the  observed  signal  for  3  June  1980 


(partition  3) ,  we  multiply  each  eigenvector  by  its  principal 
component.  The  stepwise  reconstruction  of  the  observed  vec¬ 
tor  cn  3  June  1980  data  (partition  3)  from  its  principal 
components  ar.d  eigenvectors  is  illustrated  (Fig.  15)  . 

It  is  clear  that  the  first  eigenvector  primarily 
accounts  for  variations  in  zonal  spatial  structure  vi*k 
scales  of  order  200  km  cr  greater,  and  with  significant 
amplitudes  only  beyond  80  km  offshore.  This  aspect  of  the 
scale-information  contained  in  eigenvectors  will  be  dis¬ 
cussed  in  later  sections. 

Figure  15b  illustrates  the  previously  accumulated 
signal  with  a  solid  curve,  ana  now  the  dotted  curve 
represents  the  addition  cf  the  second  eigenvector  as  modi¬ 
fied  by  the  second  principal  component.  This  mode  accounts 
for  variations  across  the  entire  transect  with  scales  of 
order  100  km.  Figure  15c  shows  the  contribution  of  the 
third  eigenvector  and  principal  component.  This  mode  has 
almost  no  effect  on  the  signal  (as  will  be  subsequently  dis¬ 
cussed,  this  happens  to  be  an  anomalous  case)  .  The  contri¬ 
butions  of  the  fourth  and  fifth  eigenvectors  and  principal 
component  are  illustrated  in  Figs.  15d  and  15e.  Comparing 
Fig.  9  with  the  final  dotted  curve  shows  that  the  data 


from  partition  three  of  VIS113  on  25  Jan  1930  ana  the  recon¬ 
structed  curve  aaree  with  the  use  of  five  eigenvalues,  and 
that  a  fair  representation  can  be  reconstructed  using  only 
the  first  two  eigenvectors. 

3.  Mean  Struct  are 

The  mean  optical  depth  (1/k(490)  =  Z90  )  transect 
profiles  for  partitions  1  through  4  are  illustrated  in  Figs. 
16  through  19  and  are  repeated  for  ease  of  comparison  in 
Figs.  24  through  27.  The  mean  vector  in  each  track 

represents  the  tendency  of  the  signal,  while  the  eigenvec¬ 
tors  scaled  by  the  principal  components  give  the  perturba¬ 
tions  of  the  mean.  In  all  four  partitions,  the  mean  value 
of  optical  depth  tends  tc  increase  with  distance  offshore. 
This  tendency  is  expected  since  the  coastal  waters  should 
contain  higher  concentrations  of  sediment  and  phytoplankton, 
especially  during  the  upwelling  season  (Traganza,  et  al. , 
1979).  There  is  a  general  lack  of  significant  eddy-like 
structure  in  the  mean  vectors  from  all  four  transects 
(although  very  lew  amplitude  perturbations  of  scale  five-to- 
ten  km  are  apparent  in  the  means). 


Distance  Offshore  km 


Reconstruction  of  Optical  Dapth  Transect  of  3 
June  1980  (Partition  3)  from  Mean  (Solid  Curve 
Panel  a)  and  Successive  Contributions  of 
Eigenvectors.  1  to  5  (Dashed  Curves)  in  Panel  a 
to  e  Respectively. 


4.  structural  Content  of  Eigenvectors  and  Principal 
Coapcnan ts  ~  ~ 

The  eigenvector  discussion  involves  many  interccm- 
parisons  of  the  partitions.  Each  partition’s  first  ten 
eigenvectors  are  plotted,  (Figs.  16  through  19  and  24  through 
27).  To  organize  the  discussion,  the  first  eigenvector  will 
be  discussed  for  all  four  partitions  before  proceeding  to 
discuss  the  second,  and  so  forth.  The  associated  principal 
components  are  also  illustrated,  (Figs.  20  through  23  and  28 
th rough  31). 

The  structure  in  the  first  eigenvector  of  each  of 
the  partitions  is  characterized  by  a  band  of  low  variability 
adjacent  to  the  coast,  and  the  structure  offshore  of  that 
band  is  dominated  by  a  scale  extending  from  there  to  the 
offshcre  end  cf  the  domain.  The  ’’node"  marking  the  onshore 
limit  of  significant  variation  in  this  mode  is  progressively 
farther  offshore,  proceeding  from  the  south  through  the  par¬ 
titions.  The  "node”  of  partition  1  begins  at  approximately 
45  km  offshore,  and  by  partition  4,  the  "node"  is  100  km 
offshore.  There  is  a  tendency  for  variance  to  decrease  in 
the  amplitude  of  the  first  eigenvector  as  the  offshore 
boundary  is  approached.  This  may  be  an  artifact  of  the 
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outer  boundary  and  should  be  investigated  further  over  a 
larger  domain  to  better  estimate  the  dominant  scale.  Parti¬ 
tion  4 ,  which  has  the  largest  spatial  extent,  shows  more  and 
larger  offshore  structure  to  beyond  180  km. 

The  associated  principal  components,  which  modify 
the  eigenvectors  before  they  are  applied  to  the  mean,  show 
the  time  variations.  Across  the  four  tracks,  the  first 
eigenvectcrs/principal  components  vary  in  phase  with  each 
other.  In  all  four  transects  there  is  a  large  difference 
between  the  first  principal  component  of  the  only  1979  image 
(early  winter)  and  those  from  1980  (early  spring).  This 
marked  difference  is  certainly  a  manifestation  of  seasonal 
variations  in  the  California  Current  system  (Pavlova,  1966  ; 
Hickey;  1979).  Host  of  the  first  principal  components  vari¬ 
ability  in  all  cases  is  otserved  in  the  1980  series  (upwell- 
ing  season)  ,  and  the  record  contains  relatively  little 
variability  in  the  1982  series  (Davidson  Current  season). 
Coherency  cf  the  variations  differs  from  oartition-to-parti- 
tion  with  nc  apparent  pattern.  The  first  eigenvector  and 
principal  component  appear  to  have  their  foundations  in  the 
offshore  seasonal  variation  and  large  scale  eddy  structure 
that  cccurs  during  the  upwelling  season.  In  the  first  mode, 


the  inshore  zone  influenced  by  upveliing  tends  to  retrain 
turbid  throughout  the  year,  whereas  the  dominant  variations 
in  optical  depth  occur  offshore  of  the  upveliing  zone. 

The  shapes  of  second  eigenvectors  from  the  four  par¬ 
titions  are  similar,  but  vary  in  an  oscillatory  fashion  from 
the  northern  partition  to  the  southern  partition.  The  first 
partition  (northern)  shovs  negative  values  beyond  approxi¬ 
mately  180  km  offshore,  and  then  small  amplitude  positive 
values  from  there  to  the  coast .  Partition  two  and  three 
depict  a  mirror  image  pattern  to  that  of  the  first  parti¬ 
tion.  Partition  four  shcvs  much  the  same  pattern  as  parti¬ 
tion  one.  The  phase  relation  in  the  principal  components 
shovs  no  pattern  between  partitions  one  and  two,  but  the 
series  for  partitions  three  and  four  both  suggest  phase 
reversal  from  the  first  partition.  This  nega tive-to-posi- 
tive-back- to-negati ve  pattern  of  behavior  weakly  suggests  a 
wave-like  meridional  oscillatory  structure,  with  an  offshore 
peak  (180  km  offshore)  in  the  vicinity  of  partitions  two  and 
three.  Resolution  of  this  meridional  characteristic  feature 
will  require  a  2-dimensional  analysis.  The  distance  of  sep¬ 
aration  of  the  partitions  suggests  a  wavelength  of  the  order 
of  120  km.  Again,  the  majority  of  the  variability  occurs 
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variability,  as  compared  tc  the  first  eigenvector. 

The  third  eigenvector  has  a  similar  behavior  for  the 
first  two  partitions.  The  perturbations  are  of  roughly  the 
same  spatial  scale  (45  km)  and  appear  to  be  in  phase.  How¬ 
ever,  the  third  and  fourth  partitions  show  an  opposite 
behavior  in  the  far  offshore  region  (beyond  183  km).  Numer¬ 
ous  smaller  scale  features  (of  the  order  of  10  km  or  less) 
are  apparent  in  this  eigenvector.  In  general,  higher  spa¬ 
tial  frequencies  become  increasingly  important  in  higher 
order  eigenvectors.  The  principal  components  show  an 
increase  in  the  variability  of  the  1982  series  with  the  wide 
range  of  variability  still  present  in  the-  1980  series.  This 
eigenvector  shows  the  largest  nearshore  amplitudes  of  all 
the  eigenvectors,  which  suggests  it  may  be  closely  linked  to 
the  nearshore  structure  of  upwelling.  There  is  little  sug¬ 
gestion  of  a  temporal  relation  evident  in  the  third  princi¬ 
pal  components  of  the  four  tracks. 

The  fourth  eigenvector  shows  an  increase  in  fre¬ 
quency  (decrease  in  wavelength)  of  the  represented  variabil¬ 
ity  scales.  Features  range  in  size  from  18  -  45  km,  with 


numerous  smaller  scale  perturbations. 
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wavelerg-h  and  range  of  variability  is  approximately  the 
same  for  all  four  partitiers.  The  principal  components  for 
partitions  one ,  three  and  four  qualitatively  suggest  coher¬ 
ency.  However,  eigenvector  behavicr  is  opposite  in  paroi- 
rion  three  and  similar  in  partition  four  when  compared  to 
partition  one.  A  prominent  feature  approximately  180  km 
offshore  ir.  the  fourth  eigenvector  of  partition  one,  seems 
to  be  shifted  outward  to  210  km  offshore  in  partition  four. 
This  time/space  relationship  between  the  structures  of  par¬ 
tition  one  and  four  again  suggests  a  meridional  oscillation 
worthy  of  future  investigation  through  a  2-dimensional 
analysis. 

The  fifth  eigenvector  shows  features  of  scale  that 
range  from  5-«2  km.  The  much  smaller  features  (less  than  5 
km)  are  not  dealt  with  as  they  are  essentially  part  of  the 
background  noise  expected  in  any  natural  system.  Little  can 
be  said  about  the  correspondence  of  the  four  partitions  with 
just  a  visual  inspection.  However,  nearshore  structure 
appears  in  the  eigenvectors  with  more  variability  than  in 
previcus  eigenvectors  of  the  1982  series.  The  1980  series 


still  demonstrates  the  largest  overall  variability,  and  the 
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large  ir.terar.nual  differences  between  the  one  scene  ir. 
and  the  1980  scenes  is  still  evident.  The  principle  compo¬ 
nents  of  partition  one  and  four  show  excellent  agreement  in 
amplitude  and  phase  for  the  first  five  time  points,  but  as 
structural  variability  decreases  with  time,  so  does  corre¬ 
spondence.  The  quantitative  correspondence  of  these  varia¬ 
tions  is  beyond  the  scope  of  this  thesis,  but  it  should  be 
investigated  in  future  analyses  of  this  data  set. 

The  sixth  through  the  tenth  eigenvectors  are  charac¬ 
terized  by  variation  of  such  high  frequency  that  little  can 
be  said  of  the  relationships  between  the  partitions.  Scales 
of  structural  features  in  these  eigenvectors  range  from  1  to 
35  km,  with  nc  suggesticn  of  a  temporal  relation  between 
partitions.  The  principal  components  show  that  the  1980 
data  again  dominates  the  variance,  but  at  these  higher  fre¬ 
quencies  the  increased  contribution  cf  the  1982  data  to  the 
total  variance  is  very  apparent.  Eecause  of  this  disorgan¬ 
ized  structure,  detailed  interpretation  is  not  attempted  for 
eigenvectors  cf  crder  greater  than  5. 
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5.  The  Joining  Of  Twc  Partitions 

Only  a  cursory  treatment  of  the  joining  function 
results  is  presented.  Partition  one  and  two  were  examined 
and  analyzed  as  per  the  joining  function  development  given 
in  Chapter  IV.  The  analysis  was  based  on  performing  the  EOF 
analysis  on  partitions  one  and  twc  combined,  and  then  com¬ 
paring  this  result  to  these  obtained  from  the  separate  EOF 
analyses  that  were  joined.  It  was  found  that  the  joining 
function  principle  components  were  within  .0001  m  cf  these 
computed  using  the  twc  partitions  as  one  data  set.  Further¬ 
more,  using  ten  degrees  of  freedom,  it  was  found  that 
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(within  2.5  X  10  )  (51) 


V  S  -  l 


(within  8.0X10  )  (52) 


This  result  is  tased  on  the  orthogonality  of  the  eigenvec¬ 
tors  as  J  is  the  matrix  cf  eigenvectors  of  the  covariance 
matrix  of  principal  components  for  the  two  partitions  com¬ 


bined. 


additionally ,  a  comparison  of  the  principal  compo¬ 


nents  yielded  (For  a  joint  sample  size  of  17) 


(A-) 


&  k  ~  h 


(within  1.4X10  )  (53) 


Represent  ing  the  principal  components  of  the  joir.ir.  g  func¬ 
tions  wicr.  Y  ,  it  fellows  that 

<✓> 

(ir)  ^  l  m  k  (within  5.0  X  10  ‘)[5«| 

Again  the  small  difference  demonstrates  the  utility  of  +  ie 
joining  process. 

The  eigenvalues  obtained  by  thus  joining  the  eigen¬ 
vectors  and  principal  components  from  partitions  1  and  2  are 
given  in  Table  III,  together  with  fractions  of  total  sample 
variance.  At  face  value,  the  first  ten  eigenvalues  account 
for  98.553  of  total  sample  variance.  Recall  however,  that 
the  input  data  were  represented  in  truncated  form,  using 
only  the  first  ten  principal  components  from  each  of  the 
partitions.  This  original  approximation  retained  only 
97.893  of  the  total  variance  computed  from  the  original 
data,  hence,  it  is  necessary  to  adjust  the  apparent  trunca¬ 
tion  cf  the  joined  result  accordingly.  The  results  of  this 
adjustment  are  given  in  the  third  column  of  Table  III  and 
show  that  assuming  that  orly  the  first  ten  eigenvectors  are 
significant  actually  leads  to  a  truncation  to  96.503  of  the 
total  sample  variance.  While  3.5%  precision  is  an  accepta¬ 
ble  level  cf  approximation  for  mest  problems  in  geophysical 


data  interpretation ,  the  effects  of  successive  *r  ur.  ca  ;  ic  r. 
must  be  giver,  careful  attention  when  applying  the  oarri- 
nicned  method  to  EOF  analysis. 

This  example  indicates  strongly  that  the  partition¬ 
ing  approach  to  EOF  analysis  can  provide  computationally 
acceptable  results  when  applied  to  satellite  image  data. 
The  example  also  emphasizes  that  proper  care  must  be  given 
to  controlling  successive  truncation  in  the  partition  join¬ 
ing  process.  It  is  left  to  future  projects  to  investigate 
questions  such  as  joining  partitions  on  the  basis  of  par¬ 
tially  intersecting  samples  to  provide  optimal  functions 
for  interpolating  satellite  data  intc  cloudy  regions,  and  to 
interpretation  of  partition  joining  functions  to  illuminate 
spatial  correlations  between  locally  important  structures 
(e.g.,  topographically  generated  mesoscale  eddies)  and  the 
dominant  structure  of  the  overall  domain  (e.g.,  that  associ¬ 
ated  with  the  evolution  of  the  synoptic  scale  upwelling 
frcnt  over  the  continental  slope  and  shelf  over  the  course 
of  the  upwelling  season) .  It  is  questions  of  this  kind  that 


address  the  ultimate  utility  of  partitioned  EOF  analysis. 
The  present  effort  is  limited  tc  preliminary  work  to  estab¬ 
lish  foundations  of  feasibility  and  procedural  constraints. 


VI.  DISCUSSION  AND  CONCLUSIONS 

Zonal  transects  of  optical  depth  (1/k(490)  m)  measured 
wirh.  the  Nimbus-7  CZCS  have  been  analysed  to  investigate 
bio-optical  structure  over  the  continental  shelf  and  slope 
off  central  California.  Samples  of  cloud  free  data  were 
selected  and  processed  for  latitudes  35-53N,  35-40N,  35-22N 

and  35-00 N .  The  data  were  observed  in  1979,  1980  ,  and  1982 
during  the  months  Hay  through  November.  The  zonal  structure 
in  these  samples  was  analysed  using  EOF's  computed  sepa¬ 
rately  for  each  section.  Meridional  variance  structure  was 
analysed  only  qualitatively  through  inspection  of  similari¬ 
ties  in  features  contained  in  EOF's  of  the  different  tran¬ 
sects  and  in  the  temporal  sequences  of  associated  principal 
components.  Finally,  the  computational  feasibility  of 
applying  partitioned  EOF  analysis  methods  to  this  type  of 
data  was  investigated  by  joining  the  EOF's  of  the  two  north¬ 
ernmost  rransects  to  form  estimates  of  the  EOF's  of  the  com¬ 
bined  spatial  domain. 

The  first  eigenvectors  for  four  zonal  transects  of  opti¬ 
cal  depth  1/(k490)  each  contained  dominant  scales  of  order 


20  0  km  or  greater,  and  accounted  fcr  between.  35  and  5  4  ~-r- 
csnt  cf  the  octal  variance.  They  are  each  also  character¬ 
ized  by  a  band  cf  low  variability  in  optical  depth  in  the 
inshore  region  influenced  by  upwelling  and  the  Davidson 
Inshore  Current.  This  band  is  confined  within  45  km  of  the 
coast  at  35  53*  N,  and  mcnoton ically  broadens  to  approxi¬ 
mately  100  km  at  35N  latitude.  This  behavior  is  possibly 
related  to  the  broadening  cf  the  continental  shelf  and  slope 
with  longshore  distance  south  of  Monterey.  Hurlburt  (1979) 
showed  that  the  topographic  beta  effect  plays  a  fundamental 
role  in  the  dynamics  associated  with  mesoscale  (order  100 
km)  longshore  variations  in  topography  by  affecting  the 
strength  cf  the  longshore  flow.  Also,  the  influence  of 
topography  can  produce  barctropic  flow  beyond  its  immediate 
vicinity.  Per  mesoscale  variations  in  coastline  geometry, 
the  coastal  currents  and  the  patterns  of  vertical  motion 
tend  to  follow  the  coastline,  but  net  with  uniform  strength. 
Coastal  current  widths  tend  to  be  narrower  than  the  scale  of 
coastline  variability.  In  these  terms,  the  meridional  vari¬ 
ation  in  scales  present  in  the  first  EOF*  s  are  qualitatively 
consistent  with  the  longshore  variations  in  bathymetry  of 
the  study  domain. 


The  second  eigenvectors  account  fcr  zonal  struccur-  win:; 
dominant  scale  cf  order  120  km ,  and  with  nearly  uniform 
amplitude  frcm  the  coast  tc  a  node  approximately  150  km  off¬ 
shore  in  partitions  1,  2,  and  4.  The  second  eigenvector  for 
partition  3  (35  20N)  is  anomalous  in  that  it  is  dominated  by 
a  zonal  waveform  with  nodes  spaced  at  approximately  80  km, 
or  roughly  half  the  dominant  scale  of  its  counterparts.  The 
reason  for  this  behavior  should  be  investigated. 

The  third  eigenvectors  are  dominated  by  scales  ranging 
from  approximately  60  to  100  km  (between  nodes) .  The  shapes 
and  scales  vary  more  strongly  frcm  partition-to-partition 
than  was  the  case  with  the  first  two  eigenvectors. 

Across  each  transact,  zonal  features  with  wavelengths 
100  km  and  greater  appear.  The  suggestion  of  an  oscillatory 
behavior  in  the  meridional  direction  needs  to  be  studied 
further.  Resolving  such  a  feature  requires  a  more  detailed 
study  involving  a  2-dimensional  analysis  cf  the  study 
domain. 

The  large  eddy  field  associated  with  the  shoreward 
boundary  of  the  California  Current  was  observed  in  the  data 
sat.  The  scales  of  this  eddy  field  were  of  the  same 
magnitude  as  the  spatial  scales  employed  in  the  partitions. 
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This  r.ac  =  ssitatsc  placing  the  partition  boundary  within  -his 
eddy  field  anc  cutting  away  seta  of  the  features.  Thn  sea¬ 
sonal  development  of  a  synoptic  scale  upwelling  front  off 
the  California  coast  is  strongly  suggested  in  the  data  and 
its  eigenvectors.  The  smaller  eddies  associated  with  this 
pattern  ranged  from  5  to  100  Jem  in  scale. 

The  convergence  of  the  eigenvalues  to  roughly  98  percent 
of  the  variance  after  the  tenth  value  was  of  particular 
interest.  This  was  true  for  all  four  partitions  and 
although  this  is  not  an  overwhelming  reduction  in  the 
degrees  of  freedem  of  the  initial  system,  it  is  significant. 

Satellite  images,  and  ether  fields  of  oceanic  and  atmos¬ 
pheric  variables,  provide  massive  data  sets.  Large  amounts 

i 

of  computer  time  must  often  be  expended  for  processing  these 
data  sets  at  even  relatively  primitive  levels.  Analyses  and 
interpretations  are,  morever,  made  difficult  by  the  sheer 
volume  of  data.  EOF  analysis  provides  a  viable  method  for 
mathematically  representing  satellite  data  fields  in  a  com¬ 
pact  and  easily  manipulated  form.  Data  transformed  using 
EOF*  s  illuminates,  and  facilitates  analysis  of,  the  time  and 
space  scales  associated  with  a  given  variable  over  the 
domain;  the  present  study  has  exercised  this  attribute  of 
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EOF's  or.  a  descriptive  level.  In  addition,  the  correct 
principal  component  representation  ci  satellite  i.r-.  j-;s  pro¬ 
vides  an  efficient  form  for  analysing  the  response  of  spa¬ 
tial  structure  in,  for  example,  optical  depth  to  forcing  by 
wind  stress  and  currents,  acting  through  a  bic-opticai 
model;  this  is  a  logical  avenue  fcr  future  research  to  build 
on  the  present  results. 

Considering  purely  computational  aspects  of  EOFs,  the 
well-known  symmetry  of  eigenvector  solutions  in  the  time  and 
space  domain  can  be  used  tc  great  advantage  in  the  analysis 
of  satellite  image  data.  The  number  of  spatial  grid  points 
in  even  the  single  trackline  partitions  of  the  present 
study  yield  large,  but  .  computationally  tractable,  scatter 
matrices.  The  larger  arrays  associated  with  2-dimensicnal 
area  partitions,  each  with  several  hundred  grid  points,  will 
clearly  exceed  sizes  admitting  direct  computation  of  spatial 
EOFs.  The  linear  algebra  and  scalings  involved  in  using  the 
smaller  time  domain  scatter  matrix  for  computation  of  space 
domain  EOF's  is  reviewed  in  Appendix  C. 

The  partitioned  method  of  EOF  analysis  illuminates  cor¬ 
relations  between  variability  in  spatially  separate  sub-re¬ 
gions.  The  present  results  demonstrate  the  computational 
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optical  depth  data.  There  is  every  reason  to  believe  -ba¬ 
the  method  may  be  equally  well  applied  to  other  CZCS  parame¬ 
ters  and  tc  infrared  imagery  of  SSI.  Further  work  in  this 
area  should  aim  to  determine  whether  the  joining  functions 
linking  EOF's  from  separate  domains  are  sufficiently  sta¬ 
tionary  tc  provide  a  basis  for  optimally  interpolating  sat¬ 
ellite  image  data  of  these  types  over  cloud-covered  areas  of 
a  particular  day's  image.  Other  applications  to  be  explored 
include  determination  of  the  extent  to  which  correlations 
between  3-dimensiona-l  in  situ  data  and  2-dimensicnai  satel¬ 
lite  data  in  small  sub-regions  may  be  extended  to  other 
parts  of  the  larger  domain  covered  by  satellite  data  alone. 


APPEND I  X  A 


SATELLITE  DATA  PROCESSING  METHODS 


A.  INTRODUCTION 

Cara  processing  was  divided  into  three  major  levels. 
Level-I  processing  includes  all  steps  required  to  take  the 
original  data  tape  to  a  Level-I  tape.  Level-II  processing 
includes  all  steps  between  a  Level-I  tape  and  a  Level-II 
tape.  Level-Ill  processing  includes  the  steps  involved  to 
take  the  Level-II  output  to  a  usable  form.  The  following 
sections  briefly  describe  the  steps  involved  in  the  three 
levels  of  processing. 

Computer  hardware  utilized  was  that  resident  at  the 
Naval  Postgraduate  School,  Monterey,  California.  The  main 
frame  computer  used  was  the  IBM  3033AP  while  the  mini-com¬ 
puter  used  was  the  Apple-* II.  Computer  software  referred  to 
in  this  section  is  either  a  system  utility  resident  to  the 
IBM  system  or  a  locally  generated  program.  Documentation  of 
the  locally  generated  programs  can  be  obtained  from: 

Dr.  J.  L.  Mueller  (Code  68My) 

Department  of  Oceanography 
Naval  Postgraduate  School 
Monterey,  California  939U3 


eft- are  involved  has 


y 


s  ystez 


veil  as  features  inserted  for 
programs  cn  other  systems  are 
mentation  carefully  prior  to 
software . 


—  d  5  d  £•  n  i  «  n  ~  f.'~  :.z  ’y  z  ~  y  ’  “ 
convieiien.ee.  Users  of  rnese 
cautioned  to  review  the  dccu- 
attempting  to  transfer  the 


3.  LEVEL -I  PROCESSING 

figure  32  is  a  schematic  diagram  illustrating  the 
processing  steps  for  Level-I  and  should  be  referred  to 
throughout  this  discussion.  The  master  tape  (raw  satellite 
data)  was  obtained  from  the  Scripps  Institution  of  Oceanog¬ 
raphy,  San  Diego,  California.  Table  IV  gives  a  summary  of 
the  master  tapes  utilized  in  this  study.  The  data  were  in 
the  ferm  of  a  standard  magnetic  tape  in  a  binary  format  with 
6250  tits  per  inch  (BPI)  .  The  tapes  were  originally  created 
using  a  Hewlett  Packard  (HE)  -  3  000  which  has  a  characteris¬ 
tic  high  order,  low  order  bit  arrangement  opposite  to  the 
IBM  system.  Therefore,  before  using  this  raw  data  in  the 
IBM  3033AP,  it  had  to  undergo  a  byte  swap  routine.  This 
byte  swap  was  accomplished  when  the  unformatted  working  bac¬ 
kup  tape  was  made  using  local  program  VISBKV.  After  the 
unformatted  backup  tape  is  made  a  variable  blocked  spanned 
(VBS)  format  tape  is  produced  using  the  system  utility 
IEBGENER. 
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Figure  32.  Level-I  Processing  Schematic  Diagram 
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con iad 


This  format  is  used  in  conjunction  with  unformatted 
statements  to  minimize  computer  time.  These  two 
tapes  serve  as  the  working  tapes  for  the  remainder  of  the 
Level-1  processing,  and  the  master  tape  is  archived. 

Using  the  VBS  formatted  tape,  a  versatec  plotter  grays¬ 
cale  is  prccuced  using  local  program  ZIPSIO.  This  program 
also  unpacks  the  Event  file  (hereafter  referred  to  as  the 
2-file)  and  writes  it  to  a  storage  disk.  The  grayscale 

depicts  the  satellite  pass  in  picture  form  for  hand  analysis 
of  landmarks.  line  numbers  and  pixel  numbers  are  taken  off 
the  grayscale  for  clear,  cloud-free  landmarks.  These  values 

are  entered  into  local  program  ZIPPIC  to  generate  a 

• 

*  P ICPRINT  * .  This  is  a  matrix  of  radiance  value  centered  on 
the  individual  landmarks  line  and  pixel  number.  These 
PICPRINTS  are  then  contoured  by  hand  (using  a  threshold 
value  of  18  counts  for  land  or  clouds)  to  determine  an  exact 
time  and  pixel  number  for  the  landmark.  The  landmark's 
latitude,  longitude,  line  number  ar.d  pixel  number  with  addi¬ 
tional  housekeeping  data  are  entered  into  local  program 


./  Jrw  r 


'■.•Wv'VlMV. 


distance  error  in  the  navigation  problem.  The  'near,  rr 
value  obtained  for  all  the  adjusted,  utilized  data  was 
approximately  1.09  n.mi..  The  final  product  of  this  step 
generates  a  set  of  navigation  parameters  that  are  used  to 
generate  a  navigation  matrix.  This  step  is  accomplished 
using  local  program  CZCSNAV2  to  generate  the  navigation 
matrix  and  NAVDUMP  to  write  the  navigation  matrix  (here¬ 
after  referred  to  as  the  G-file)  to  a  temporary  formatted 
tape. 

The  E-file  is  copied  from  its  temporary  disk  storage  to 
the  temporary  tape  as  the  G-file.  Additionally  the  Data 
file  (hereafter  referred  to  as  the  L-file)  is  first  reversed 
from  its  tcttcm-to-top  orientation  to  a  top- to- bottom  orien¬ 
tation  using  local  program  BACKWARD.  This  program  also  puts 
the  L-file  to  the  previously  mentioned  temporary  storage 
tape.  Finally,  these  files  on  the  temporary  tape  are  copied 
to  a  Level-I  tape  using  the  system  utility  IEBGENER.  The 
only  difference  between  the  temporary  tape  and  the  final 
Level-I  tape  is  that  the  L-file  is  copied  into  an  unformat¬ 
ted  file  which  will  aid  in  the  speed  of  further  processing. 
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C.  LEVEL -II  EEOCESSI  SG 

figure  33  is  a  schematic  diagram  illustrating  the 
processing  steps  for  Level-II.  The  Lavei-I  taps  generated 
by  the  steps  discussed  in  the  previous  section  is  the  input 
tape  for  this  processing.  Only  the  L-file  is  affected  by 
the  Level-II  processing  as  the  E-file  and  G-file  are  copied 
straight  to  the  Level-II  tape  using  the  system  utility 
IEBGENES.  The  L-file  is  used  to  generate  output  for  calcu¬ 
lating  the  proper  values  of  the  Angstrom  coefficient  for 
each  scene.  This  is  done  using  local  program  CZPARMS2  and 
an  assumed  value  for  the  Angstrom  coefficient.  Chapter  III 
Section  C.2.  discusses  the  importance  and  method  of  finding 
these  values.  Next,  the  computed  Angstrom  coefficients  with 
the  L-file  are  rerun  through  CZPARJJS2  to  regenerate  the 
L-file.  This  regeneration  involves  taking  the  raw  counts  of 
each  channel  and  applying  the  bio-optic  algorithms  discussed 
in  Chapter  III  to  produce  values  for  chlorophyll  and  K.  At 
this  point  the  adjustment  algorithm  discussed  in  Chapter  V 
has  no t  been  applied.  Analyses  cf  the  initial  Level-II  out¬ 
put  precipitated  the  need  for  the  corrective  algorithm, 
which  was  then  applied  during  Level-Ill  processing. 
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Figure  33.  Level-II  Processing  Schematic  Diagram 
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LEVEL- III  PROCESSING 


Figure  3*J  illustrates  the  Level-Ill  processing  steps. 
The  Level-Ill  processing  basically  takes  the  data  obtained 
in  Level-II  and  marries  it  to  the  navigation  matrix  gener¬ 


ated  during  Level-I  processing.  Using  the  four  designated 
tracks  (Fig.  5  ),  the  coastal  starting  points  from  each 
track  were  entered  into  local  program  TED DUMP  to  provide  the 


navigation  blcck  of  the  track  origin.  The  G-file  contains 
data  (in  latitude  and  lcngitude  values)  every  sixteenth 


pixel  and  sixteenth  line.  Once  the  origin  block  is  estab¬ 


lished  the  exact  line  and  pixel  was  interpolated  using  local 
program  FINDPIX  on  the  Apple  II.  With  this  starting  point 
local  program  DATA4  was  entered  to  generate  every  1  km  along 
each  track  an  associated  line  and  pixel  number  which  was 


then  converted  into  the  appropriate  data  values.  This  out¬ 
put  was  written  to  storage  for  later  processing.  It  was 
hare  that  the  adjustment  algorithm  was  applied,  producing 
the  final  version  of  the  data  in  a  navigated  form. 
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Figure  34.  Level-Ill  Processing  Schematic  Diagram 
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Figures  36  through  39  show  the  plots  generated  by  par- 
plot  and  the  partitioning  given  in  Table  IV.  The  objective 
of  the  partitioning  was  to  find  the  most  complete  data  over 
time  and  space  possible  given  ten  samples. 

Cata  frcm  these  four  partitions  was  then  entered  into 
local  program  PARDAT.  This  program  applied  most  of  the  con¬ 
ditioning  to  the  data  set.  Only  the  K(4  9Q)  data  was  uti¬ 
lized  from  this  point  on  although  this  program  could  be 
easily  altered  to  focus  on  another  optical  parameter.  The 
data  were  searched  to  find  good  points  and  bad  points  and  a 
control  arrangement  for  later  use  was  made.  The  raw  data 
were  scaled  and  inverted  to  produce  1/K(490)  values  in 
meters.  Next  the  data  were  averaged  by  every  fourth  point 
to  smooth  out  noise  features.  At  this  point  data  strings 
with  gaps  existed  for  each  applicable  scene.  Next  a  linear 
interpolating  routine  was  applied  to  obtain  continuous  data 
at  each  time  point.  Finally,  the  data  for  all  scenes  and 


tracks  were  combined  and  written  to  disk  in  a  single  data 
file.  To  this  pcint  the  conditioning  applied  has  consisted 
of  partitioning  the  data  into  four  partitions,  rejecting 
incomplete  scenes,  scaling  the  K  (490)  values,  averaging  the 
data  by  every  four  values,  and  applying  linear  interpolation 
to  fill  in  the  remaining  gaps. 

The  conditioned  data  was  then  plotted  using  local  pro¬ 
gram  PARPLT  and  the  DISSELA  system  utilities.  The  plots 
generated  are  figures  7  through  10  and  were  used  in  Chapter 
V  Section  B  to  discuss  the  data  and  its  relationship  to  the 
regicral  oceanography. 

The  final  steps  of  the  data  conditioning  involved  appli¬ 
cations  of  the  EOF  analysis  techniques.  Local  program  PEOF1 
produced  eigenvalues,  eigenvectors,  and  principle  components 
for  each  partition  and  plotted  the  output.  Figures  11 
through  14  and  16  through  31  are  the  plots  produced.  The 
eigenvalues,  eigenvectors  and  principle  components  were  all 
written  to  disk  for  later  use.  The  local  program  JE0F2  was 
designed  to  produce  the  joining  function  that  related  parti¬ 
tion  one  to  partition  two. 
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Figure  38.  Partitioning  Scheme  for  Track  Three  (35  22  K) 
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Figure  39.  Partitioning  Scheme  fcr  Track  Four  (35  53  N) 
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The  desired  EOF’s  are  those  in  the  space  domain ,  which 
for  satellite  dara,  is  dimensioned  much  larger  than  the  time 
dimension.  It  is  possible  to  significantly  expedite  compu¬ 
tations  by  computing  eigenvalues,  eigenvectors  and  principal 
components  using  the  smaller  covariance  matrix  of  the  time 
domain,  and  to  then  scale  these  results  to  obtain  the  eigen¬ 
vectors  and  principal  components  in  the  space  domain.  The 
algebraic  basis  for  this  approach  is  reviewed  in  this 
appendix . 


Consider  the  followinc  convention  for  dimension  nota¬ 


tions 


Spac< 


n  •  1  , .  . . ,  N 


Time  .  .  .  .  m  ~  1 , .  . . ,  U 

EO F  (order)  .  .  .  .  k  —  1 , . .  . , K 

where  K  <  min(M,N).  As  before  the  raw  data  matrix  is  given 


KJ 


of  M  rows  by  N  columns.  The  sample  mean  in  the  space  domain 


is  given  by 


'  Y  -  -  ~ •  '  • 1  •  jn 


f  -  M 


M 


u-=  1 


mn 


(56) 


of  dimension  1  X  N.  The  centered  data  matrix  in  the  space 
domaiE  is  given  by 


l 


■[«■'] 


(57) 


From  Chapter  IV  the  sample  space  ccvariance  matrix  (of  siz 
N  X  N)  is 


M-l  F  F 


(53) 


and  the  sample  time  covariance  matrix  (centered  in  space  and 
of  size  a  X  M)  is  now  defined  as 


T 

a. 


M-l  F  F  • 


(59) 


If  A  and  e  are  an  associated  eigenvalue  and  eigenvector 

K.  “  K 

respectively  cf  S,  then 


S  e,  =  A,  e,  ^  0  , 

^  -k  k  -k 


(60) 


since 


A^  #  0  and  i  0 


(61) 


Then 


F  e,  f  0  , 

'v-  -k  ’ 


(62) 


and  therefore 


M-l  t 


t  r 


l 


A  .  f  , 

k  ^  -  k. 


(*■*) 


(noting  the  commutative  property  of  and  ).  Thus,Ak  is 
an  eigenvalue  of  T  with  associated  eigenvector  (in  the  time 
domain) 


^k 


l  2i 


(6U) 


where  amk  is  the  k'th  principal  component  (in  the  space 
domain)  at  time  point  m. 

We  now  consider  the  eigenvalues  x  and  eigenvectors  £, 

k.  —  K 


of  T,  normalized  such  that 


M 


6ij 


m  =  1 

with  time  domain  principal  components 


(65) 


Sk  =  l  Sk  5  k  -  1, . • . ,K  > 


(66) 


where  the  •  significance  order  K  is  selected  either  to  retain 
an  arbitrary  fraction  of  the  total  sample  variance,  or 
according  to  one  of  the  more  objective  statistical  selection 
rules  discussed  by  Preisendorfer ,  et  al.  (1981). 

The  principal  components  in  the  time  and  space  domain 
have  the  property  that 


N 


n  *  1 


M 


CnJ 


2Z  a«i 


(M-l)  Xi  5 


(67) 


m  =  1 
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5-=call  that  F  was  centered  on  f  in  the  space  domain,  ar. 
that  both  S  ar.d  T  were  therefore  normalized  by  the  facto 

-A 

(1/M-1)  . 


Ey  analogy  tc  equation  (63) 

mTT  I1  I  ST  Sk  *  xk  f  Ik  '  xk  Ik 


where  5.  is  now  seen  to  be  an  eigenvector  of  S,  but  o 

**  K.  ^ 

length 


(M  -  1 


lk  ’ 


rather  than  being  an  orthcnormal  eigenvector  of  length  1 
The  two  are  thus  related  by 


Sk  = 


(principal  components  in  the  time  domain  scale  to  orthonor 
mal  eigenvectors  in  the  space  domain)  , 


and  by  similar  arguments 


(M  -  l)\ 


k  ^k 


(orthonormal  eigenvectors  in  the  time  domain  scale  to  prin 
cipal  components  in  the  space  domain)  . 
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